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ABSTRACT 


Radar  imaging  is  an  area  of  tremendous  interest  as  radar-based  systems  are 
perhaps  the  only  all-weather,  long  range  remote  sensing  systems.  However,  radar’s 
continued  utility  and  application  in  wide-ranging  areas  is  fundamentally  dependent  on  the 
ability  to  produce  high-quality,  artifact- free  imagery.  To  date,  the  use  of  radar  to  identify 
and  image  moving  objects  remains  of  great  interest,  and  it  is  well  known  that  motion  in 
the  scene  gives  rise  to  mispositioning  or  streaking  when  target  motion  is  not  addressed. 
Many  techniques  have  been  developed  to  handle  moving  objects,  but  these  techniques 
make  use  of  the  start-stop  approximation,  in  which  a  target  in  motion  is  assumed  to  be 
momentarily  stationary  while  it  is  being  interrogated  by  a  radar  pulse. 

A  new  linearized  imaging  theory  that  combines  the  spatial,  temporal  and  spectral 
aspects  of  scattered  waves  has  been  developed.  This  thesis  studies  the  performance  of 
this  new  imaging  scheme  via  analysis  of  the  point  spread  function.  It  is  shown  that  the 
imaging  PSF  localizes,  and  is  translation  invariant  in  phase-space.  It  is  also  shown  that 
the  behavior  of  the  imaging  system  is  dependent  on  the  aperture  geometry  and  choice  of 
radar  waveforms  transmitted. 
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I.  INTRODUCTION 


A.  RADAR  IMAGING 

Sensor  systems  that  can  detect,  locate  and  identify  targets  at  great  distances  and  in 
all  kinds  of  weather  have  well-recognized  utility.  Radar-based  systems  are  capable  of 
performing  such  functions  over  long-range,  all-weather  and  all-day  (day  or  night)  unlike 
many  optical  counterparts.  This  is  because  the  wavelength  of  radar  signals  makes  them 
relatively  unaffected  by  atmospheric  and  weather-induced  attenuation.  Herein  also  lies 
the  problem:  image  resolution  depends  upon  signal  wavelength.  The  resultant  side  effect 
is  that  resolution  suffers  and  radar-based  imaging  will  be  fundamentally  less  effective 
compared  to  other  systems  based  upon  relatively  shorter  signal  wavelengths. 
Nevertheless,  radar  waves  scatter  strongly  from  objects  whose  size  is  on  the  same  order 
as  the  wavelength.  This  means  that  radar  is  sensitive  to  objects  whose  length  scales  range 
from  centimeters  to  meters.  Many  objects  of  interest  are  in  this  range. 

Many  research  and  development  efforts  have  been  expended  to  improve  the  use  of 
radar  for  imaging.  It  has  been  reported  that  resolution  on  the  order  of  centimeters  can  be 
achieved  with  spotlight-mode  operation  [15].  While  the  developed  techniques  have  been 
instrumental  in  making  radar  imaging  a  more  viable  application,  the  parallel  advent  of 
fast  digital  computers  and  leaps  in  computational  prowess,  have  enabled  practical 
systems  to  be  fielded.  This  has  generated  great  activity  in  areas  such  as  interferometric 
synthetic  aperture  radar  (SAR)  and  bistatic  SAR,  due  to  interest  in  operational  concepts. 
The  key  to  success  in  all  of  these  applications  is  the  ability  to  produce  high-quality, 
artifact-free  imagery. 

B.  IMAGE  ARTIFACTS 

The  subject  of  radar  (and  radar  imaging)  has  evolved  considerably  since  the 
1950s.  Consequently,  many  implicit  assumptions  are  embedded  within  current  imaging 
techniques.  These  simplifications  or  system  limitations  are  of  interest  as  they  are 
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frequently  also  responsible  for  the  image  artifacts  discussed  in  the  literature.  The  start- 
stop  approximation  is  central  to  the  thesis  and  holds  the  fundamental  difference  between 
the  new  imaging  scheme  developed  in  [2]  and  current  techniques. 

1.  Scattering  Model  -  Born  Approximation 

Historically,  scattering  models  have  been  built  upon  a  foundation  which  views 
targets  as  being  composed  of  simple,  non-interacting  point  scatterers.  When  the  imaging 
scene  consists  of  multiple  targets,  accounting  for  multiple  scattering  (in  which  waves 
bounce  back  and  forth  between  the  targets)  can  become  extremely  complex.  For  this 
reason,  the  ‘weak  scatterer’  or  Bom  approximation  is  typically  invoked.  A  useful 
consequence  is  that  the  imaging  problem  becomes  linear  and  more  tractable. 

However,  the  Born  approximation  may  result  in  image  artifacts  as  actual 
“multiple-scattered”  waves  are  received  later  than  expected  from  the  single-scattering 
model.  These  waves  are  then  interpreted  as  having  come  from  a  more  distant  target 
structure.  In  addition,  increased  capabilities  in  radar  resolution  and  sensitivity  and 
parallel  advancements  in  computational  technology  have  exposed  the  inadequacies  of 
such  simple  scattering  models.  With  the  ability  of  radar  systems  to  resolve  and  image  the 
components  (subscatterers)  that  lie  within  the  support  of  a  traditional  target,  it  is 
inevitable  that  new  radar  target  models  that  can  describe  these  multiple  scattering  effects 
accurately  (whether  multiple  targets  or  multiple  components  of  a  single  target)  are 
required.  Current  approaches  typically  utilize  correction  techniques  to  enhance  imaging 
algorithms  based  on  the  low  order  ‘weak  scatterer’  model. 

2.  Measurement  Systems 

There  are  inherent  limitations  in  any  data  that  are  measured  by  radar  systems. 
Such  data  are  usually  corrupted  by  noise  and  are  always  limited  by  bandwidth. 
Consequently,  inducing  unwanted  image  artifacts  and  reduced  image  resolution. 
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3. 


“Start-Stop”  Approximation 


The  radar  uncertainty  principle  bounds  the  accuracy  with  which  time  delay  and 
frequency  shift  can  be  simultaneously  estimated  using  a  single  pulse.  This  limitation  may 
be  mitigated  by  the  use  of  follow-on  pulses  to  further  probe  target  characteristics.  Most 
modem  radars  use  a  train  of  high  range-resolution  waveforms  to  enable  pulse-to-pulse 
velocity  estimation.  Such  waveforms  allow  the  use  of  the  start-stop  approximation, 
which  assumes  that  the  target  is  stationary  during  the  time  the  radar  pulse  illuminates  the 
target,  and  moves  only  between  pulses.  When  individual  pulses  are  short,  and  targets  of 
interest  are  small  and  move  slowly  compared  to  the  speed  of  light,  this  is  almost  always  a 
good  approximation.  Under  these  circumstances,  good  velocity  estimates  are  obtained 
from  the  range-rate  (two  accurate  range  measurements,  separated  by  a  time  interval). 
Unfortunately,  no  Doppler  information  is  actually  measured  by  the  systems.  Range  and 
range-rate  information  are  then  used  to  form  spatial  images. 

In  order  to  achieve  the  necessary  signal-to-noise  ratio  (by  coherent  integration), 
long  pulse  trains  are  required.  Consequently,  during  the  coherent  processing  interval,  the 
stationary  target  assumption  may  no  longer  hold.  Scatterers  may  move  to  different 
positions,  and  when  they  move  beyond  a  resolution  cell,  the  result  is  a  blurred  image. 

C.  MOTIVATION  -  IMAGING  MOVING  TARGETS 

In  recent  years,  a  number  of  attempts  to  develop  imaging  techniques  that  can 
handle  moving  objects  have  been  proposed.  Space-time  adaptive  processing  (STAP)  is  a 
signal  processing  technique  that  was  originally  developed  for  detecting  slow-moving 
targets  using  airborne  radars.  STAP  uses  multiple-element  antenna  array,  together  with 
real-aperture  imaging  techniques  to  produce  ground  moving  target  indicator  (GMTI) 
images  [23],  This  method  described  in  [24]  combines  SAR  (designed  to  image  stationary 
scenes)  and  GMTI  processing  to  develop  new  techniques  for  detecting  slow-moving 
surface  targets  that  exhibit  start-stop  maneuvers.  The  method  described  in  [11]  discusses 
the  extension  to  non-sideways  looking  array  radar. 
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The  method  described  in  [10]  describes  the  concept  of  velocity  synthetic  aperture 
radar  (VSAR).  VSAR  is  a  multi-element  SAR  system  involving  conventional  processing 
to  form  an  image  at  each  element;  except  the  image  phases  are  preserved  and  compared 
to  deduce  target  velocity.  However,  VSAR  processing  assumes  that  the  scatterers  remain 
in  a  single  resolution  cell  throughout  the  integration  time. 

Time-frequency  methods  for  signal  analysis  allow  decomposition  of  the 
frequency  spectrum  of  time -varying  signals  in  shorter  time-frames,  and  provide  an  added 
dimension  for  examining  the  dynamic  behavior  of  the  signal  as  it  varies  over  time.  The 
size  and  shape  of  the  weighting  window  function  can  be  altered  to  fit  the  specific  analysis 
requirement  for  the  signal.  In  recent  years,  further  extension  of  this  time-frequency 
method  to  ISAR  has  been  advocated  in  the  area  of  image  generation,  motion 
compensation  and  micro-Doppler  target  vibration  studies  [Refs.  8,  20  and  21].  Chen  [20] 
has  reportedly  used  the  time-frequency  transfonn  (TFT)  to  construct  inverse  synthetic 
aperture  radar  (ISAR)  images  of  targets  with  high  rotation  rate  by  helping  to  overcome 
the  migration  of  individual  scattering  points  from  one  range  cell  to  another. 

However,  all  of  these  techniques  rely  on  making  the  approximation  that  a  target  in 
motion  is  assumed  to  be  momentarily  stationary  while  interrogated  by  a  radar  pulse.  It  is 
apparent  that  the  utility  of  an  imaging  approach  that  can  accommodate  target  motion 
during  the  imaging  process.  This  constitutes  the  focus  of  the  thesis. 

D.  OBJECTIVE 

The  objective  of  this  thesis  is  to  study  a  linearized  imaging  theory  developed  by 
Cheney,  Borden  [2]  for  an  imaging  region  containing  moving  objects.  The  physics 
involved,  and  the  approach  for  addressing  image  artifacts  common  to  targets  moving  in 
an  unknown  fashion  will  be  discussed.  Analysis  of  the  corresponding  point-spread 
function  will  be  carried  out  to  examine  the  localization  in  phase-space  (position-velocity). 
The  behavior  (position,  and  velocity  resolution)  of  the  imaging  system  will  also  be 
examined  under  different  aperture  geometries  and  transmitted  waveforms. 
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E.  THESIS  ORGANIZATION 

The  thesis  is  organized  in  the  following  manner: 

Chapter  II  serves  to  develop  the  appropriate  radar  imaging  theory  by  looking  at 
the  basic  building  blocks  necessary  for  imaging  -  radar  operating  principle,  scattering 
models,  correlation  reception  -  and  how  these  components  translate  to  imagery.  Finally, 
the  thesis  includes  an  overview  of  the  inverse  nature  of  the  radar  imaging  problem  to 
provide  better  appreciation  of  issues  pertaining  to  practical  measurement  systems  and 
how  information  is  extracted  from  measurements  to  obtain  the  best  estimates  of  the  target 
object. 

Chapter  III  applies  the  concepts  put  forth  in  Chapter  II  to  the  imaging  scenario. 
The  discussion  here  is  based  on  the  work  of  Cheney,  Borden  [2]  in  imaging  moving 
targets  from  scattered  waves.  The  effect  of  multiple  moving  targets  on  the  scattering 
model  is  described.  The  approximations  invoked,  specifically  slow-moving,  narrow-band 
and  far- Held,  will  be  examined.  A  scattering  model  to  be  used  in  the  simulation  program 
will  be  presented. 

Chapter  IV  discusses  the  radar  ambiguity  function,  and  the  relationship  between 
estimation  errors  and  the  transmitted  waveform.  Particular  emphasis  will  be  given  to 
image  quality  analysis  and  the  derivation  of  the  point-spread  function  to  quantify  the 
performance  of  the  imaging  scheme  developed  in  Chapter  III. 

Chapter  V  examines  further  implications  of  the  transmitted  wavefonn  through  the 
pulse  compression  perspective.  In  particular,  the  Chirp  signal  is  discussed. 

To  demonstrate  the  behavior  of  the  imaging  scheme,  a  simple  simulation  program 
is  developed  for  modeling  the  point-spread  function.  Chapter  VI  discusses  the  results  of 
the  simulation,  specifically  the  effect  of  different  aperture  geometries  and  transmitted 
waveform  characteristics.  It  also  concludes  and  summarizes  the  major  findings  of  the 
thesis  work  and  provides  recommendations  for  future  work. 
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II.  IMAGING  THEORY 


Prior  to  deriving  the  radar  data  model  for  imaging  moving  targets,  Chapter  II 
examines  the  formation  of  spatial  images.  It  begins  by  examining  the  physics  behind 
extracting  range  and  range-rate  (Doppler  shift)  information  by  a  radar  system.  A  simple 
range-Doppler  imaging  example  is  presented.  The  chapter  concludes  with  a  brief  look  at 
the  inverse  nature  of  the  imaging  problem  and  characteristics  of  measurement  systems, 
and  from  there  some  of  the  issues/challenges  pertaining  to  radar  imaging. 

A.  BASIC  RADAR  PRINCIPLES 

Most  fielded  radar  systems  are  pulsed  systems,  as  opposed  to  stepped-frequency 
systems.  Radar  systems  measure  the  strength  of  the  backscattered  field  and  the  round  trip 
delay  time  of  transmitted  signal  pulse  reflected  from  distant  objects.  Since  the  radar 
pulse  travels  at  the  speed  of  light,  it  is  practical  to  use  the  measured  time  delay  for  the 
round  trip  to  determine  the  range  to  the  reflecting  object.  This  is  the  basic  principle  of 
radar  theory  and  we  can  detennine  the  range  R  from  R  =  ct  /  2  (where  t  is  the  round  trip 
time  delay). 

Several  concepts  will  be  useful  in  appreciating  how  radar  signals  convey 
information  for  the  imaging  output.  The  first  is  concerned  with  the  kind  of  signals  that  a 
radar  uses  to  generate  the  waveforms  it  transmits.  Radar  signals  are  real-valued  functions 
of  time  but  it  is  convenient  to  express  the  wavefonns  as  complex-valued  functions. 
Radar  systems  are  able  to  measure  these  signals  (in  both  amplitude  and  phase)  as 
functions  of  time.  The  imaginary  part  of  the  complex  signal  provides  the  phase  of  the 
waveform  revealing  considerable  information  about  the  target. 

Radar  waves  are  electromagnetic  waves  governed  by  Maxwell’s  equations.  The 
waveforms  generated  from  radar  signals  satisfy  the  wave  equation  and  the  radiation 
condition.  These  properties  ‘code’  the  echo  (transmitted  wave  scattered  from  a  radar 
target)  waveform  with  important  target  specific  information  such  as  range,  speed  and 
bearing. 
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The  resolution  in  the  down-range  direction  is  governed  by  the  pulse  width,  and 
correspondingly  its  pulse  bandwidth.  Typical  carrier  frequency  operation  bands  occur 
within  the  low  attenuation  atmospheric  windows  ~  1-40  GHz.  The  frequency  band  of  the 
system,  in  practice,  signifies  that  radar  signals  are  always  of  finite  duration.  Intuitively,  a 
higher  bandwidth  will  imply  finer  resolution  for  imagery  in  the  range  direction.  More 
about  the  concept  of  “resolution”  will  be  addressed  in  Chapter  IV. 

B.  SCATTERING  OF  ELECTROMAGNETIC  WAVES 

Understanding  the  behavior  of  electromagnetic  waves,  specifically  propagation 
and  scattering  of  these  waves,  is  essential  to  appreciating  the  infonnation  that  can  be 
extracted  from  the  transmitted  signals.  The  next  two  sections  will  elaborate  this  further. 

The  discussion  of  [18]  describes  a  simple  one-dimensional  scattering  model  for  a 
moving  PEC1  plate  in  free  space  (a  good  approximation  for  typical  radar  wave 
propagation).  An  elementary  consequence  of  Maxwell’s  equations  is  that  an 
electromagnetic  wave  (in  the  far-field)  must  have  the  electric  field  magnetic  field 

A 

and  wave  vector  (direction  of  wave  propagation)  k  mutually  perpendicular.  For  the 
transmitted  wave  with  field  components  orthogonal  to  the  x  direction  (i.e.  in  the  y-z 
plane),  the  boundary  condition  on  the  PEC  results  in  the  following  relationship  between 
scattered  and  transmitted  electric  field: 

fr,,,gu+m)ic)=-»cnt-w)ic) 

where  the  scattered  field  consists  of  a  left-travelling  wave  (in  the  negative  x  direction). 
The  position  of  the  moving  plate  is  given  by  R(t) .  Assuming  for  simplicity  that  the  plate 
is  undergoing  linear  motion2  (with  velocity  v  in  the  x-direction),  so  that 


1  A  perfect  electrical  conductor  (PEC)  allows  the  charges  to  move  freely  and  instantaneously  in 
response  to  a  field;  consequently  the  fields  inside  a  PEC  are  zero.  The  resulting  boundary  conditions  are: 
the  tangential  components  of  the  electric  field  must  be  zero;  and  the  tangential  components  of  the  magnetic 
field  are  related  to  currents  flowing  on  the  surface  of  the  PEC. 

2  Even  if  the  motion  is  not  linear,  the  linear  expression  R(t)  =  x+vt  can  be  considered  to  be  the  first  two 
terms  in  a  Taylor  expansion  that  is  valid  for  small  t  (i.e.,  for  short  pulses). 
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R(t )  =  x  +  vt 


x=R( at  boundary ) 


and  making  use  of  the  substitution  u  =  t  +  R(t )  /  c ,  we  can 


solve  for  t  in  terms  of  u  (at  the  boundary): 


u  -  R  /  c 

t  = - 

1  +  v/c 


and  obtain  the  expression  for  the  scattered  field 

+  R(t)  /  c)  =  +  xl  c-RI  c)-RI  c] 


(2.1) 


where  a  is  the  Doppler  scale  factor 


1-v/c 
1  +  v/c 


(2.2) 


For  a  transmitted  field  consisting  of  a  signal  waveform  s( t)  mixed  with  a  carrier 
wave  cos (co  t)  ,  the  scattered  signal  received  at  the  antenna  (x  =  0)  is 

prec(t)  *  s(t  -  2  R  /  c)cos((U0  [(1  -  2v  /  c)(t  -R/  c)-R/  c])  (2.3) 

We  have  used  the  fact  that  s(t )  is  slowly  varying  to  approximate  the  argument 
a(t  -  R  /  c)~  R  /  c  by  t  -2R  /  c  ,  and  also  the  fact  that  for  a  scatterer  moving  slowly,  so 
that  vie  is  small.  In  addition,  we  can  expand  the  denominator  of  Equation  (2.2)  in  a 
geometric  series  and  approximate  a  by  1-2 vie.  The  multiplier  of  t  in  the  argument  of 
the  cosine  function,  is  the  carrier  frequency  (relative  to  the  transmitted  carrier  frequency) 
shifted  by  an  amount 

2v 

v  = - c o0  (2.4) 

c 


This  is  the  Doppler  shift.  As  in  the  stationary  scatterer  case,  the  received  signal  is  time- 
delayed,  but  here  it  is  also  a  Doppler-shifted  version  of  the  transmitted  signal. 


C.  CORRELATION  RECEPTION  -  RADAR  DATA  MODEL 


Detection  of  signals  in  noise  is  fundamental  to  radar  operations.  A  radar  detects 


the  presence  of  an  echo  signal  reflected  from  a  target  and  extracts  information  about  the 
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target  for  uses  such  as  imaging.  The  conservation  of  energy  requires  that  the  wave 
intensity  from  a  finite-sized  antenna  decreases  in  strength  as  the  inverse  square  of  the 
distance.  Consequently,  the  round  trip  power  of  the  signal  received  by  the  same  antenna 
for  small  targets  decays  by  a  factor  of  1  /  R4  relative  to  the  signal  that  was  transmitted. 
For  typical  radar  operational  distances,  the  received  power  can  be  as  low  as  10'  Watts 
consequently.  Thermal  noise  in  the  receiver  presents  a  problem  for  target  detection,  and 
a  greater  problem  for  image  formation. 

While  increasing  the  power  of  the  transmitter  (or  pulse  integration)  can  help 
increase  the  signal-to-noise  ratio  (SNR),  there  are  practical  limits  to  relying  solely  on 
these  approaches.  Key  amongst  these  limits  is  the  decrease  in  maximum  unambiguous 
range,  and  also  the  fact  that  targets  move  during  data  collection  for  pulse  integration, 
altering  the  phase  of  the  scattered  field.  It  should  be  pointed  out  that  the  ensuing 
discussion  is  premised  on  the  correlation  comparison  being  performed  coherently  such 
that  the  phase  of  the  transmission  signal  is  preserved  in  the  reference  signal. 

In  practice,  a  solution  involves  transmitting  long  coded  pulses,  together  with 
appropriate  signal  processing  techniques  called  matched  filtering  or  correlation 
reception.  This  is  known  as  pulse  compression  (see  Chapter  V).  The  correlation  receiver 
(or  matched  filter)  is  an  important  example  of  a  radar  signal  processor  for  detection  of 
desired  signals.  The  output  of  the  matched  filter  is  the  cross-correlation  function  of  the 
received  signal  and  the  expected  received  signal  (which  is  related  to  the  transmitted 
signal).  Hence,  the  mathematical  equivalence  between  the  correlation  receiver  and  the 
matched  filter  receiver,  and  it  is  possible  to  implement  the  matched  filter  as  a  correlation 
process. 

A  radar  obtains  information  about  a  target  by  comparing  the  received  echo  signal 
with  the  signal  that  was  transmitted.  Correlation  reception  presents  an  intuitive  approach 
to  develop  the  standard  radar  data  model  used  in  imaging.  The  signal  processing  problem 
is  that  of  optimal  detection  in  additive  noise.  Assuming  Gaussian  noise  statistics, 

s  (t') 

maximum  likelihood  processing  (conditioned  on  the  measured  signal  rec  )  selects  the 
best  guess  estimate  for  the  radar  signal  associated  with  the  field  reflected  from  a  target 
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g  ( 

scan v  ’  that  maximizes  the  conditional  probability.  Reference  [19]  shows  that  such 


processing  (averaged  over  all  time)  reduces  to  finding  the  function 


,(0 


that  minimizes 


J  kciO\2  dt'+  J  \sscjn[  dV-  2Re  J  Srec(t'yscatt(t')dt' 


s  (t')  s  (t') 

The  objective  is  to  determine  scatt  ’  from  the  random  measurements  of  rec v  . 

The  first  two  terms  are  system-scenario  specific.  Searching  for  the  best  minimizing 

function  hinges  on  the  last  term  and  can  be  made  more  efficient  by  restricting  the  search 

to  a  few  parameters  derived  from  some  model.  In  this  case,  the  natural  model  is  based  on 

s  ( t') 

the  scattering  interaction  between  the  interrogating  field  and  the  target.  If  tnc v  ’  denotes 
an  incident  pulse  transmitted  by  the  radar,  then  the  linear  radar  (weak-scatterer)  scattering 
model  follows  by  superposition 


SXJ0  =  jj  -  T)em~T)d rdv  (2.6) 


where  p{o,  r)  is  the  target  reflectivity  density  for  a  radar  system  that  measures  two 
parameters  (range  and  radial  velocity).  This  reflectivity  density  is  defined  in  such  a  way 
that  p(o,T)dzdo  is  proportional  to  the  field  reflected  from  the  target  at  range  between 
cr/2  and  c(T  +  dr)l  2  with  Doppler  shift  between  v  and  o+do. 

From  Equation  (2.6),  the  scattered  field  is  naturally  modeled  by  the  two 
parameters  r  and  u  .  The  simplified  search  space  will  be  parameterized  by  the  time  shift 
r  and  the  frequency  shift  u .  Considering  the  last  term  of  Equation  (2.5),  and 
maximizing  for  all  r ,  o  in  the  scattered  model,  the  correlation  integral  is 

oo 

r/(u, t)  =  |  srec{t')s]nc(t'-  v )e~iv(t'~T)dt'  (2.7) 


The  correlation  integral  seeks  out  the  component  of  received  signal  that  “matches”  the 
time-delayed,  frequency-shifted  version  of  the  transmitted  signal.  Unlike  signals  will  be 
suppressed  consequently;  correlation  receivers  attempt  to  find  the  r  and  v  that 
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maximize  the  real  part  of  Equation  (2.7).  The  expected  output  of  the  correlation  receivers 
is  obtained  by  substituting  Equation  (2.6)  into  the  correlation  integral  of  Equation  (2.7) 


rj(o,r)  =  JJJ  p(v',r')s.nc(t'-  T)s*nc(t'-  r')eiv(t'  T)e^wi,'~T,)dt' dr' do'  (2.8) 


Recall  that  the  signal-processing  problem  was  detection  of  signals  in  noise,  hence 
the  received  data  sr  (t)  consists  of  both  the  scattered  signal  and  a  random  noise 

component.  The  correlation  noise  term,  which  is  the  correlation  integral  between  random 
noise  and  the  time-delayed,  frequency-shifted  version  of  the  transmitted  signal,  has  been 
dropped  from  Equation  (2.8).  The  correlation  noise  term  which  measures  how  well  the 
scattered  signal  matched  the  random  noise  will  usually  be  small  for  coherent  systems. 
Rearranging  Equation  (2.8), 


oo 

7(u,r)  =  ||  p(o',t')z(v-v',t - 

-00 


i— (u+u')(t-t') 

r')e2  dr'du' 


(2.9) 


where 


CO  ,  . 

Z(o,  t)  =  J  s.nc(t'-  -  T)s*nc(t'+  -  z)eivt  dt'  (2.10) 

-oo 

Equation  (2.9)  is  the  standard  radar  data  model  and  expresses  the  output  of  the  correlation 
receiver  as  the  convolution  of  p  and  z  (UP  to  a  phase  factor).  The  function  z(v>  r) 
defined  by  Equation  (2.10)  is  the  radar  ambiguity  function.  Radar  imaging  is  closely 
related  to  the  radar  ambiguity  function  that  characterizes  the  accuracy  that  the  target 
position  and  radial  velocity  are  estimated  from  the  radar  data  (see  Chapter  IV). 

D.  ONE-DIMENSIONAL  (HIGH  RANGE  RESOLUTION)  IMAGING 

One-dimensional  imaging  involves  the  generation  of  range  profiles,  which  can  be 

used  to  articulate  target  substructure.  For  high  range  resolution  (HRR)  techniques,  the 

transmitted  pulse’s  instantaneous  range  resolution  is  smaller  than  that  of  the  target,  and  as 

it  sweeps  across  the  target  it  sequentially  excites  the  target’s  scattering  subelements 

which  re -radiate  energy  back  to  the  receiver.  When  these  scattering  subelements  are  non- 
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interacting  and  point-like,  the  scattered  pulse  will  be  a  sum  of  damped  and  blurred 
images  of  the  incident  pulse,  which  are  shifted  by  time  delays  that  are  proportional  to  the 
subelement’s  range  [1].  Figure  1  shows  an  example  of  one-dimensional  images  created 
by  HRR  radar  systems.  It  can  be  observed  that  the  profile  displays  additional  peaks 
outside  of  the  target  support,  clear  evidence  of  the  problem  associated  with  a  single,  non¬ 
interacting  scatterer  approximation. 


y  (meters) 

0  10  20  30  40 


Figure  1 .  Example  of  a  range  profile  from  a  B-727  jetliner.  The  top  view  (with 
orientation  at  time  of  measurement)  is  displayed  beneath  [1], 
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The  signal  received  from  a  distribution  of  non-interacting,  stationary  targets  can 
be  modeled  as 

oo 

s,.Jt)  =  J  p(  x ')sjt  -  z')dz '+  n(t )  (2.11) 

-oo 

where  n(t)  represents  an  additive  and  random  noise  term.  Application  of  a  matched 
filter  to  obtain  the  radar  data  model  used  for  imaging 

oo  oo  oo 

r/(z)  =  J  srec(t')s*nc(t'~  z)dt'  =  |  s*iic(t'~  t)  J  p(z')s.nft'-  z')dz'dt'  +  noise  tenn 

—oo  —oo  —oo 

oo 

=  JJ  s*nc(t'~  z)s.nc(t'~  z')dt' p{z')dz'  +  noise  tenn 

—00 

Making  the  substitution  t"  =  t'—z',  and  neglecting  the  noise  term  (negligible  with 
matched  filtering)  gives 

oo 

7(0=  J  jO-rXryr'  (2.12) 

-oo 

where  x  is  the  autocorrelation 

oo  oo 

X(f)=  J  s*(t"—  z)s(t")dt"  =  J  s*(t')s(t'+  z)dt'  (2.13) 

-oo  -oo 

Equation  (2.12)  represents  a  one-dimensional  image  and  shows  how  the  true  distribution 
p  is  related  to  the  radar  data  model.  The  image  is  the  convolution  of  p  with  x,  which  in 
high-range  resolution  (HRR)  imaging  is  also  commonly  referred  to  as  the  point-spread 
function  (PSF).  The  point  spread  function  obtains  its  name  from  the  fact  that  if  p  consists 

of  a  single  point  S(t),  then  the  image  is  7(0  =  J  xif  ~  t')8(t')dt'  =  x(0  ■  Thus,  x 

quantifies  the  degree  to  which  a  single  point  appears  spread  out  in  the  HRR  image. 

The  down-range  profile  can  be  affected  by  a  variety  of  factors  such  as  target 
aspect  angle,  position  of  the  scatterers  or  masking  of  scatterers  by  other  parts  of  the 
target.  Additionally,  while  the  use  of  short  pulses  enhances  the  resolution,  it  concurrently 
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leads  to  large  bandwidth  requirements.  Wide  bandwidth  can  increase  system  complexity 
and  increase  the  likelihood  of  interference  from  other  emitters  in  the  electromagnetic 
spectrum.  A  short-pulse  waveform  also  provides  less  accurate  radial  velocity 
measurement,  a  natural  consequence  from  the  properties  of  the  radar  ambiguity  function 
(more  details  on  radar  resolution  and  the  ambiguity  function  in  Chapter  III).  An 
important  implementation  limitation  to  highlight  for  practical  radar  applications  is  the 
required  high  peak  power  to  transmit  short  pulses  over  long  ranges.  High  peak  power 
transmission  can  be  problematic,  especially  at  high  frequencies  (due  to  the  small 
waveguide  dimensions  and  the  small  area  of  the  anodes  of  the  microwave  sources). 

Radar  target  recognition  using  only  range  profiles  has  limited  applications.  This 
is  because  a  range  profile  will  not  be  able  to  distinguish  cross-range  target  structures.  All 
scatterers  located  at  the  same  distance  from  the  radar  will  reflect  energy  back  to  the  radar 
with  the  same  time  delay.  Hence,  when  the  radar  illuminates  many  distinct  targets  at  any 
instant,  meaningful  interpretation  of  on-scene  target  (multiple)  disposition  cannot  be 
derived  based  on  a  set  (single)  of  range-only  data. 

E.  TWO-DIMENSIONAL  IMAGING 

For  more  effective  target  interpretation,  additional  infonnation  is  needed  on  top  of 
the  range  profile.  This  information  can  be  in  the  form  of  high-resolution  cross-range 
profile,  Doppler  profile,  or  simply  the  “triangulation”  of  range  profiles. 

As  a  simple  illustration  for  extension  of  the  radar  imaging  concept  to  two 
dimensions,  “triangulation”  of  different  sets  of  range  profiles  will  be  examined  (Figures  2 
and  3).  This  approach  enables  the  determination  of  cross-range  target  structure  while 
using  only  HRR  radar  systems  and  relies  on  collecting  multiple  sets  of  range  profiles 
from  different  target  orientations,  processing  them,  and  subsequently  synthesizing  a  two- 
dimensional  image. 

Consider  a  set  of  three  point  targets  with  the  radar  located  at  the  same  distance 
from  targets  2  and  3  (see  Figure  2).  When  the  geometry  is  oriented  as  shown,  the  return 
echo  will  only  indicate  two  targets.  Hence,  ambiguity  exists  when  targets  lie  along  the 
bands  of  constant  range  from  the  radar. 
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Target  1 

mam 


Bands  of  constant 
range  from  radar 


Target  3 


o 


o 


Target  2 


Single  pulse  ambiguity 


Figure  2.  Ambiguous  scenario  from  a  single  radar  pulse  [19]. 


With  multiple  sets  of  data  from  different  directions,  triangulation  will  allow  for 
gradual  buildup  of  the  relative  positions  of  the  three  targets  (Figure  3).  The  range 
profiles  are  swept  in  the  cross-range  direction  to  form  bands  of  constant  range  from  the 
radar,  which  represent  the  possible  locations  of  the  target  scatterers.  These  bands  are 
then  superimposed  and  the  crossing  points  are  used  to  determine  the  scattering  center 
locations.  The  success  of  such  a  scheme  hinges  fundamentally  on  the  ability  to  correlate 
the  various  constant  range  bands.  Otherwise,  the  intersection  of  the  swept  lines  will  be 
erroneous,  resulting  in  image  artifacts. 


Multiple  pulse  “imaging” 


Target  1 

o 


Target  3 


o 

Target  2 


o 


Figure  3.  Cross-range  information  obtained  from  range  profiles  [19]. 
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The  signal  due  to  a  moving  target  is  a  time-delayed,  Doppler-shifted  version  of 
the  transmitted  signal.  The  measured  echo  signal  from  a  single  point-like  moving  target 
can  be  expressed  as 


s(t)  =  ps,-(t  -  z)e~iv(t  T)  +  n(t ) 


(2.14) 


With  two  unknown  parameters  to  estimate,  namely  the  delay  z  and  the  Doppler 
shift  v ,  a  single  matched  filter  will  no  longer  suffice.  A  set  of  matched  filters,  one  for 
every  possible  Doppler  shift,  will  have  to  be  applied 

oo  oo 

J  lift  -  t')sr  (t')dl'  =  pj  hft  —  t')sinc(t'— z)e  r]dt'  +  noise  tenn 

—oo  —oo 

To  estimate  z  and  u ,  find  their  values  where  1 77 1  takes  on  its  maximum.  For  a 
distribution  of  non-interacting  moving  targets,  the  output  from  the  filter  bank  is 

00 

r/(z,o)=  J  s*nc(t'~  z)e~w(,'~T)srec{t')dt' 

-00 

00  00 

=  J  s*nc(t ’—  z)e~w{t  ~r)  ||  p(z',u')sjnc(t'-  z')e~w(,~T)dz'do'dt' 

-00  -00 

00 

=  JJJ  s*nc(t'-  z)e-iv(t'-T)sJt'-  z'y-^'-^dV p(z' ,v')dz' do'  (2.15) 

-00 

00 

=  JJJ  4e(t"+  z'-  zy-W'^'-^s^Oe^dt " p(z\v')dz'dv' 

-00 

00 

=  JJ  z\v-  u')e~w(T~T"> p(z\u')dz' do' 

-00 

where 


Z(r,u)=  J  s*nc(t"+  z)e  w(n s(t")dt"  =  J  s*nc(t'+  z)sinc(t')e  w,'dt'  (2.16) 

-00  -00 

is  the  radar  ambiguity  function.  Similar  to  the  one-dimensional  case,  the  ambiguity 
function  can  be  interpreted  as  an  imaging  point-spread  function  for  range-Doppler 
imaging.  Analysis  of  the  ambiguity  function  determines  the  fidelity  of  the  range-Doppler 
imaging  process,  and  consequently  the  accuracy  of  target  range  and  velocity  estimation. 
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1.  Range-Doppler  Spatial  Image 


With  estimates  for  range  and  range-rate  from  radar  returns,  this  information  can 
be  interpreted  to  form  spatial  images.  The  following  treatment  illustrates  this,  provided 
the  relative  motion  between  the  antenna  and  the  target  is  known.  Consider  a  simplified 
two-dimensional  (2D)  geometry  between  the  radar  and  target  as  shown  in  (Figure  4). 
The  target  is  assumed  to  have  translational  and  rotational  motion  only  in  the  two- 
dimensional  plane  relative  to  a  stationary  radar  platfonn. 

Consider  a  rotating  target  [18].  For  simplification,  assume  that  any  overall 
translational  motion  is  removed. 


Figure  4.  The  geometry  for  locating  a  point  target  known  to  be  rotating  [18]. 

Consider  a  point  target  rotating  counter-clockwise  in  a  geometry  in  which  the 
antenna  is  in  the  plane  of  rotation  (Figure  4).  The  coordinate  origin  coincides  with  the 
center  of  rotation;  and  the  center  of  rotation  is  a  large  distance  R()  from  the  radar  antenna. 
The  large  range  approximation  means  that  the  incident  wave  is  approximately  a  plane 
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wave,  consequently,  all  horizontal  points  are  (approximately)  at  the  same  range.  If  the 
target  is  positioned  at  (x,y)  =  a(cos#,sin#) ,  then  its  range  is 

R  =  RQ  +  y=RQ  +  a  sin#  (2.17) 

If  its  angular  velocity  is  $=  d6 /  dt ,  then  its  velocity  vector  is  sin#,  cos#),  its 
tangential  speed  is  v  |=  a&,  and  the  down  range  (y)  component  of  the  velocity  is 
$=|  v  |cos#  =  aStosO.  However,  since  x  =  a cos#,  this  down-range  velocity  can  be 

written  $=  x§*,  where  x  is  the  cross-range  component  of  the  position.  Note  that  for 
x>  0 ,  the  down-range  velocity  is  negative,  whereas  for  x  <  0 ,  it  is  positive. 

This  means  that  since  a  radar  system  can  measure  range  (Ru  +  v)  and  range-rate 
(&=x$),  it  can  determine  the  (x,y)  location  of  the  point  in  space  (provided  #  is 
known).  In  particular,  the  radar  measures  the  time  delay  t  =  2R  /  c  and  the  Doppler  shift 

°D  = - ^0= - 

c  c 

which  means  that  the  coordinates  of  the  point  are  given  by 

{x’y)={-vtb  (218) 

This  shows  how  the  range-rate  radar  measurements  can  be  used  to  estimate 
Doppler  and  cross-range  information.  Further,  the  simplification  introduced  highlights  to 
the  difficulties  in  practical  implementation  of  the  radar  scattering  model.  Estimates  of  the 
scattered  field  requires  any  variations  in  range  R  to  be  estimable.  While  it  is  reasonable 
to  assume  (in  the  one-dimensional  imaging  methods)  R  as  constant  over  the  data  set,  the 

two-dimensional  imaging  schemes  invariably  will  not.  R  will  generally  vary  due  to  the 
target  (translational)  motion  that  occurs  during  measurements  made  while  the  target 
rotates.  The  effect  of  failing  to  correctly  account  for  this  translational  motion  will  be  the 
introduction  of  phase  errors  into  these  data,  and  a  source  of  image  corruption. 
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F. 


RADAR  IMAGING  -  AN  INVERSE  PROBLEM 


A  discussion  of  imaging  theory  will  not  be  complete  without  addressing  the 
inverse  nature  of  the  imaging  problem  and  its  associated  issues.  As  the  thesis  is  focused 
primarily  on  developing  and  analyzing  the  point  spread  function  rather  than  imaging,  the 
intent  of  this  section  is  to  highlight  the  key  characteristics  of  such  inverse  problems  - 
issues  pertaining  to  practical  measurement  systems  and  how  infonnation  is  typically 
extracted  from  measurements  to  obtain  the  best  estimates  of  the  target  object. 

The  radar  can  be  well-approximated  by  a  linear,  shift-invariant  system.  This  is  an 
important  characteristic,  as  it  states  that  the  object  function  f(x)  at  each  point  x  has  a 
corresponding  output  image  m(x ) .  Further,  the  system  will  measure  objects  the  same  no 
matter  when  they  start  (or,  equivalently,  where  they  are  located).  Radar  imaging  will  not 
be  meaningful  if  the  delayed  version  of  the  target  looks  different.  A  general  formulation 
of  the  imaging  problem  can  be  achieved  by  considering  an  imaging  system  represented 
by  the  functional  operator  (“kernel”)  k  that  describes  how  the  measurement  system 
works 

m  =  Kf  (2.19) 

The  direct  problem  in  radar  imaging  refers  to  the  mapping  from  the  target  to  the 
quantities  that  can  be  measured  by  the  radar.  Consequently,  the  inverse  problem  is 
concerned  with  ‘reproducing’  the  original  target  from  given  data  and  knowledge  of  the 
direct  problem  (in  the  case  of  radar  imaging,  it  will  be  radar  measurements  and 
knowledge  of  the  scattering  model). 

1.  Well-Posed  and  Ill-Posed  Problems 

One  of  the  issues  for  inverse  problems  is  their  ill-posed  nature.  An  ill-posed 
problem  is  one  whose  solution  is  not  unique  and/or  does  not  exist  for  any  data  and/or 
does  not  depend  continuously  on  the  data.  In  the  case  of  a  band-limited  system,  the 
solution  of  the  inverse  problem  is  not  unique.  This  is  because  the  imaging  system  does 
not  transmit  information  about  the  target  at  frequencies  outside  the  band  of  the  measuring 
system. 
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Consider  discretized  measurement  systems  in  which  the  measurements  are 
collected  at  specific  points,  and  so  Equation  (2.19)  can  be  written  in  matrix  form.  And 
the  basic  goal  of  getting  a  good  estimate  /  from  the  measurements  m  can  be  notionally 
represented  as 

/  =  K'm  (2.20) 

It  can  be  observed  that  k  ,  which  represents  the  measurement  system,  imposes  certain 
properties  on  the  measurement  system.  In  general,  the  matrix  k  is  N  x  M  (  M  >  N  )  as 
the  object  function  is  infinitely  dimensional  and  will  not  be  confined  to  the  limitations  of 
the  measurement  system.  Consequently,  Equation  (2.20)  has  more  unknowns  than  there 
are  linearly  independent  equations  and  so  the  system  cannot  have  a  unique  solution. 

Since  the  object  space  (M-dimensional)  is  typically  larger  than  the  measurement 
space  (N-dimensional),  there  exists  what  is  known  as  the  nullspace  which  consists  of  all 
vectors  f-icf  or,  physically,  all  the  things  that  cannot  be  measured.  Nullspace 
accounts  for  measurement  artifacts  and  is  detennined  by  the  kernel  k  . 

The  accepted  approach  for  solving  inverse  problems  which  are  ill-posed  is  to 
search  for  approximate  solutions  satisfying  additional  constraints  based  on  the  physics  of 
the  problem  [22],  This  set  of  approximate  solutions  corresponding  to  the  same  data 
function  is  the  set  of  objects  with  images  close  to  the  measured  one  and  is  expressed  in 
what  are  known  as  the  normal  equations 

KTKf  =  KTm  (2.21) 


2.  Data  Reconstruction  -  Regularization 

Equation  (2.21)  specifies  the  set  of  conditions  that  must  be  satisfied  by  the  least- 
squares  solution  /  which  can  be  rewritten  as 

/  =  (kt  kYxkt  m  (2.22) 

Recall  that  the  measured  data  suffers  from  finite  dimension  and  noise  contamination. 
These  are  problems  since  {kt k)  is  poorly  conditioned  and  the  effect  of  noise  is 
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magnified  (i.e.  errors  increase)  as  the  dimensionality  of  the  problem  gets  larger.  Two 
techniques  will  be  highlighted  here,  namely  the  singular  value  decomposition  (SVD)  and 
Tikhonov  regularization,  to  mitigate  the  effects  of  noise  when  estimating  / . 

From  linear  algebra  methods,  the  matrix  {ktk)  xkt  can  be  expressed  in  terms  of 
its  SVD  [19]: 

(ktk)~1kt  =  UDVt  (2.23) 

where  the  matrices  U,  V  contain  the  eigenvectors  of  (ktk),  and  D  will  have  diagonal 

elements  of  the  fonn  1  /  (where  A.  are  the  corresponding,  non-zero  eigenvalues).  A 

consequence  of  the  kernel  k  being  bounded  is  that  the  eigenvalues  necessarily  form  a 
sequence  that  gets  arbitrarily  close  to  zero  as  the  number  of  measurements  increase.  The 
diagonal  elements  of  D  get  large  when  the  higher  dimensional  data  are  used,  and  if  the 
data  are  contaminated  by  noise,  from  Equations  (2.22  and  2.23)  that  noise  will  also  be 
multiplied  by  a  large  value.  Hence,  one  simple  approach  to  mitigate  the  noise  effect  is  to 
truncate  the  SVD  representation  when  A.  falls  below  some  chosen  threshold  value.  The 

corresponding  diagonal  element  is  set  to  zero  and  this  method  is  known  as  the  truncation 
filter. 

Tikhonov  regularization  modifies  ( KTKfxKT  to  a  better  conditioned  form  by 
adding  a  fixed  scalar  a  to  the  original  nonnal  equations.  Equation  (2.22)  can  then  be 
rewritten  in  the  form 

/  =  {k1 k  +  al)~l  kt m  (2.24) 

When  a  =  0,  the  original  nonnal  equations  are  returned,  so  it  can  be  expected  that  small 
values  of  a  will  not  change  the  system  description  too  much. 

After  introducing  a  ,  the  diagonal  elements  of  D  in  Equation  (2.23)  become 

£>„  =  -&-  (2.25) 

A.  +  oc 

i 
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When  A  ?  a  ,  D  «  1  /  J~A.  (as  before).  When  /L.  =  a ,  however,  D  «  /  a  and  will 

not  become  arbitrarily  large.  Consequently,  the  noise  problem  is  controlled. 
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III.  IMAGING  DATA  MODEL 


This  chapter  develops  an  approach  to  imaging  moving  objects  from  scattered 
waves.  Many  of  the  ideas  and  concepts  introduced  in  the  preceding  chapter  will  be 
extended  to  the  three-dimensional  case  in  order  to  develop  the  imaging  model. 

A.  INTEGRAL  EQUATION  APPROACH  TO  SCATTERING 

We  have  seen  in  the  preceding  chapter  that  for  the  one-dimensional  scattering 
problem  from  a  moving  plate,  the  signal  that  returns  to  the  radar  is  a  time-delayed  and 
frequency-shifted  version  of  the  transmitted  signal.  This  is  also  the  case  for  three- 
dimensional  scattering  from  a  moving,  point-like  target.  A  point  target  is  the  traditional 
model  used  for  developing  imaging  analysis  and  will  also  be  used  for  the  thesis.  The 
representation  of  a  target  as  points  is  appropriate  because  for  high-frequency  waveforms 
scattered  from  smooth  conducting  surfaces,  the  net  scattering  effects  from  the  surface  is 
equivalent  to  the  contribution  from  the  entire  series  of  points  known  as  specular  points 
[19]. 

Applying  traditional  electromagnetic  approaches  to  scattering  will  typically 
involve  solving  a  differential  wave  equation  and  imposition  of  appropriate  boundary 
conditions  at  the  surface  of  the  scatterer.  In  practice,  this  process  is  difficult  to 
implement  for  there  is  no  coordinate  system  that  can  represent  the  boundary  conditions 
for  complex  targets  we  encounter  in  practical  scattering  applications. 

Hence,  it  is  convenient  to  refonnulate  the  scattering  problem.  It  turns  out  that 
there  exists  an  equivalent  integral  equation  for  every  differential  equation,  and  this 
integral  equation  can  intrinsically  incorporate  any  boundary  conditions  in  the  solution. 

We  take  the  scattering  body  to  be  located  in  free  space  and  describable  by  an 
index  of  refraction  n(x) .  Consider  a  monochromatic  wave,  with  hannonic  time 
dependence  (any  signal  can  be  built  from  these  individual  frequency  component  via 
superposition).  The  reduced  wave  equation  for  the  field  (Helmholtz  equation)  is 
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V2x¥(x)  +  k2n\x)x¥(x)  =  0  (3.1) 

where  k  is  the  wave  number  in  free  space  (note  that  n(x )  ^  1  on  the  scatterer).  Splitting 
the  field  VF  into  the  incident  field  vF.m  representing  the  incoming  wave,  and  the 


scattered  field  *¥  scaU  and  simplifying,  Equation  (3.1)  yields 

+  ^..Jx)  =  kHl  -  n\x))(v,Jx)  +  V„„«) 

=  l  'P . (vi) 

where  the  scattering  density  p(x)  =  k2  -  n2 (x)^  is  a  source  factor  completely 

detennined  by  the  index  of  refraction  of  the  scatterer  and  is  zero  outside  of  the  scatterer. 

A  solution  to  Equation  (3.2)  can  be  obtained  based  on  the  application  of  Green 
functions,  and  we  can  write  the  solution  as 

+  (3.3) 

This  is  a  Lippmann-Schwinger  equation.  It  can  be  observed  that  the  problem  of  solving 
equation  (1.1)  has  been  reduced  to  finding  the  Green  function  Gk(x\x)  in  free  space  (k 
denotes  a  specific  frequency).  The  scattered  field  is  also  written  as  an  integral  over  the 
scatterer  D  since  p(x)  =  0  otherwise,  and  hence  automatically  satisfies  the  required 
boundary  conditions  for  a  known  target. 

It  is  appropriate  to  make  several  additional  observations  with  respect  to  the 
Lippmann-Schwinger  equation: 

(1)  There  are  generally  many  functions  that  can  serve  as  Green  functions 
Gk(x',x)  in  free  space.  While  there  are  no  bounding  surfaces  (other  than  those  defined 
by  the  scattering  object),  we  are  interested  in  Gk(x\x)  that  obey  the  Sommerfeld 

radiation  condition  lim Gk{x',x)  oc  - 
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(2)  The  inverse  scattering  problem  -  i.e.  that  of  determining  p  given  *P  and 
¥ ,  can  also  be  approached  using  the  Lippmann-Schwinger  equation  as  a  model  for 
'F  . 

scatt 


B.  LINEARIZED  DATA  MODEL  (TIME- VARYING  SYSTEMS) 


Enforcing  the  Bom  (“weak  scatterer”)  approximation,  where  the  object 
distribution  (points)  is  not  self-interacting,  multiple  scattering  events  are  treated  to  be 
negligible  in  comparison  with  the  primary  scattering  events.  The  resulting  linearized 
(approximate)  data  model  is 

(*>><** ')'*'.«  (*  W  (3.4) 

This  solution  is  one  for  the  Helmholtz  equation.  For  time-varying  systems,  the 
time-dependent  wave  equation  and  the  time-domain  Green  function  will  have  to  be 
incorporated  within  the  model 

^  ,  ,  .  ,  ,  ,  .  8(t-t'-\x'-x\lc)  _ 

Gk(x  ,x)  -»  g(x  ,x;t  ,t)  = - - — — - - -  (3.5) 

k  4k  \  x  -  x  | 

For  a  signal  s  (t)  transmitted  from  position  y  at  starting  time  —T  ,  the  incident  field  at  x 
is 


y/inSx^= 


sy(t  +  Ty-  \x-y\/c) 
4n\x-y\ 


(3.6) 


and  so  the  scattered  field  y/  (y,z,t)  at  time  t  and  position  z  obeys  (Born 
approximation)  is 


r(v,z,0  =  JJ 


cc8(t-t'- 

z  —  x' 

1  lc)%(f+T- 

\*'-y\ 

'J  4k 

z  —  x'\ 

4k 

X 

'~y 

p(x')d3x'dt'  (3.7) 


where  denotes  the  second  time  derivative  of  sv(t)  and  enters  as  a  consequence  of  the 
second  time  derivative  in  the  wave  equation. 
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c. 


REFLECTIVITY  FUNCTION  FOR  MOVING  TARGETS 


The  reflectivity  function  p(x)  or  scattering  density  first  introduced  in  Equation 
(3.2)  is  completely  specified  by  the  target  scatterer.  It  is  a  scale  factor  of  the  received 
signal  strength  and  an  important  characteristic  of  radar  imaging.  For  a  moving  target,  a 
time-varying  reflectivity  function  p(x,t)  will  be  used.  For  an  imaging  scene  containing 
multiple  moving  targets  (or  scatterers),  the  scattering  model  can  be  modified  as  follows: 
Let  p  (x  -  vt)d3xd3v  be  the  scatterers  in  the  volume  element  cVxd  'v  of  phase  space 
centered  at  position  x  and  velocity  v.  Choose  T  to  be  such  that  the  scattering  density 
pv  =  pv(x)  describes  those  targets  moving  with  velocity  v  at  time  t=0  (i.e.,  the 
transmitting  antenna  is  activated  at  time  t  =  —T  ).  The  spatial  scatterer  density  centered 
at  time  t  and  position  x  is 

P(x)  =  \pv(x-vt)d3v  (3.8) 

and  the  scattered  field  obeys 


|  z-x' 

1  /c) 

JJJ  4^| 

z  —  x'  \ 

k(t'+T-\x'-y\/Q 
An\x'-y\ 


(3.9) 


p  (x'~  vt')d3vd3  x' dt 


Making  the  change  of  variables  x'  —>  x  =  x'—vt'  (i.e.,  change  to  a  frame  of  reference  in 
which  the  scatterer  p  is  fixed)  to  obtain 


S(t  -t'  —  |  x  +  vt'-  z\/c ) 
An  I  x  +  vt'-  z  I 


'+T -\x  +  vt'-y\/c)  ,  , 

—  p  (x)d  vd  xdt' 


(3.10) 


An  |  x  +  vt'-  y  \ 


The  physical  interpretation  is  that  the  wave  that  emanates  from  y  at  time  —T  encounters 
a  target  at  time  t' .  This  target,  during  the  interval  [0,C],  has  moved  from  x  to  x  +  vt' . 
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The  wave  scatters  with  strength  pv(x)  and  then  propagates  from  position  x  +  vt'  to  z, 
arriving  at  time  t. 

Simplifying  the  notation,  and  writing 

R(t)  =  x  +  vt-z,  R  =|  R  |,  R  =  R/R 

then 


V, 


.(**)= JJJ 


Sjt-t'-  Rx  z(t')/  c) 

4^(0 

B(t'+T  -R  ( t')/c ) 

x  ^ -  ( x)d\dixdt ’ 

4  n  R  (O 


(3-11) 


The  variable  t’  in  the  argument  of  the  delta  function  appears  in  several  terms  making  it 
difficult  to  carry  out  the  integral  over  t’.  To  address  this,  let  /'  =  t  t  (t)  denote  the 

implicit  solution  of 


t-t'-  RrXt')/  c  =  0,  i.e.,  t-txv(t)-Rxz(tx,v(t))/c  =  0 


Then,  performing  the  integration  over  V  yields 

&(txv(t)  +  Ty-Ry(txv(t))/c) 

V scanty =  j|  ,.\7„  —  )-  7^rPM)d  vdx 


Substituting  t ,  (t)  =  t  -  R  Xt -  ,(0)  /  c  for  the  first  occurrence  of  t  (t)  in  this  last 


equation  allows  for  a  more  symmetric  version 


Vs, 


,(v,z,t)  =  JJ 


S \(t  +  T -R  (t  (i t))/c-R  ( t  {t))/c\ 

y  \  y  x,z v  x,v v  77  x,y v  x,v v  77  / 


p  (x)d3vd3x  (3.12) 


D.  FURTHER  SIMPLIFYING  APPROXIMATIONS 

Equation  (3.12)  applies  to  very  general  circumstances.  It  applies  to  rapidly 
moving  targets  (where  the  start-stop  approximation  is  not  valid),  and  is  appropriate  for 
arbitrary  transmitted  waveforms. 
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However,  Equation  (3.12)  is  also  a  manifestation  of  retarded-time  problems, 
which  are  invariably  complex.  The  complications  arise  as  a  consequence  of  the 
definition  of  retarded  time  txv(t),  which  is  obtained  implicitly  from  the  relation 

t  (t)  =  t-  R  (t  (tj)  /  c.  A  more  tractable  result  can  be  obtained  when  the  scatterers 

are  known  to  be  moving  in  such  a  way  that  R  (t)  can  be  expanded  in  a  Taylor  series 

around  t=0  and  approximated  by  retaining  only  the  terms  linear  in  t  -  i.e.,  when  the 
scatterer  is  “slow  moving”. 


1.  Slow-Mover  Approximation 


Assume  that  (|  v  1 1)  and  (|  v\2  t2  x  com  /  c)  are  much  less  than  |  x  -  z  |  and 
|  x  -  y  |  (where  eo  denotes  the  maximum  angular  frequency  of  the  transmitted  signal 
5  ).  In  this  case, 

y'  ’ 


RxJt)  =\z-(x-  vt)  |=  Rxz(  0)  +  iix  z(0)g?t  + ... 

where  R  ( 0)  =  jc-z,  i?v  (0)  =|  7?  (0)  | ,  and  Rx  z(0)  =  Rx  .(0)/  Rx  ,(0).  Substituting 
this  result  into  the  definition  of  retarded  time  yields 

U0^-(y  0)  +  RX0)9>Tu(t))/c 

„  <-yo)/c 

~  i  +  iix  z(0)9>/c 


Inserting  this  approximation  into  the  result  for  y/  ( y,z,t )  yields 


¥ 

'  si 


t 


y  V  x’v  L 


c  - R .  (0  )/c  +  T 


x,y 


1 


pi(x)d3vdix  (3.13) 


(4?r)2  Rx  z(0)Rx  y(0) 


where 


\-Rxy(0)g>/c 
l  +  f?v.(0  )g;/c 


(3.14) 
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is  the  Doppler  scale  factor.  Similar  to  the  traditional  Doppler  shift,  it  is  closely  related  to 
the  velocity  component  in  target-transmitter  direction. 


2,  Slow-Mover  and  Narrow-Band  Approximation 


Typical  systems  are  narrow-band  for  which  the  transmitted  signal  is  of  the  form 


s(t)  =  §Mt)e 


-ico  t 


where  §/<(t)  is  slowly  varying  (as  a  function  of  /)  in  comparison  with  e  Wy  (where  co  is 
the  carrier  frequency  for  the  transmitter  at  position  y).  In  this  case,  the  time  derivatives  of 
s  are  dominated  by  the  e  w>y  factor,  and  St(t)  «  -co2§j(t)e  '°>f' .  Moreover,  since  §j(t) 
is  slowly  varying,  further  approximations  can  be  made 


§/a (ax  J)e 


■  §/c {t)e 


Equation  (3.13),  under  the  slow-mover  and  narrow-band  then  becomes 
Yscatl{y,z,t)  =  -\\- 


j  irn  -iO)av  t 

co  e  x,ve  y  ’ 
y 


(4xYRMRy(0) 


%(t  +  T:  -  (Rx  _(0)  +  R  v(0))  /  c^y>v(x)d3vd3x 


where 


(3.15) 


(3.16) 


R  (0 )-cT  +  a  R  (0) 

x,y  v  7  y  x, v  x,z  v  ' 


/  c 


(3.17) 


3.  Slow-Mover,  Narrow-Band  and  Far-Field  Approximation 

When  the  transmitter-to-target  and  target-to-receiver  distances  are  large  in 
comparison  with  the  scene  dimensions,  then  Ijc  +  vtl  and  I  x  +  vt  I2  xco  tc  can  be 

1  ~  l  i  i  i  max 

assumed  to  be  much  less  than  either  |  z  |  or  |  y  \ .  And  the  expansion 

Rx  _(0  =|  z-(x-  vt )  |=|  z  |  -£o(x  +  vt)  +  ... 
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can  be  applied  (and  similarly  for  Rx  (t)).  Substituting  the  expansions  into  previous 
approximations  (Equation  (3.16))  then  yields 


■■ 


®  ff  i(p  -ico  a  t 

- —?■ -  e^e 

(4^)2|z||y|JJ 

x  ,%(t  +  Ty  -  (|  z  |  -igv+ 1  y  |  -dgv)  /  c  y  (x)d3vd3x 


(3.18) 


where  now 


(O  =  co 

~  X.V  V 


\y\-y&-cT+a{\z\-tgc) 


/  c 


(3.19) 


Note  that  in  Equation  (3.16),  the  expansions  for  R  (/) ,  R  (/)  are  used  differently.  In 
the  amplitude  factor  the  approximation  R  _(t)  =|  z  -  (x  +  vt)  |«|  z  |  holds  since 

|  z  |?  |  igi  x  +  vt)  | ;  in  the  phase  factor,  however,  the  product  (£gv  x  comm  /  c)  can  still  be  a 

large  fraction  of  2 n ,  and  must  be  retained.  In  addition,  the  signal  envelope  §/  can  still 
vary  significantly  over  |  z  |  -iigv  despite  being  slowly  varying,  hence  the  first  order  term 
is  retained. 

The  Doppler  scale  factor  also  reduces  to 

au~\  ^~7~  ~  1  ~  (P+  ®)9>  / c  (3-20) 

1  +  £C y/c 

where  the  binomial  expansion  is  used  since  |  v  |  tc  =  1.  The  quantity  co  x  (]  is  the 

Doppler  shift  (where  /?  =  -(p+  ®)cy  t  c )  and  is  observed  to  be  dependent  on  the 
“bistatic”  vector  ( y>+  i5). 

Inserting  T  =\y\/c  and  setting  k  =  coy  / c 

<PXtV  =  K  |  z  |  -ky(y>+  fi)^x  +  (rg[z  -  x))v  /  c]  (3.21) 
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E. 


IMAGING  VIA  A  FILTERED  ADJOINT 


In  developing  the  correlation  receiver,  it  was  observed  that  the  “fit”  between  two 
(complex-valued)  functions  is  given  by  the  cross  correlation.  For  correlation  reception, 
the  correlation  integral  seeks  out  the  component  of  received  signal  that  “matches”  the 
time-delayed,  frequency-shifted  version  of  the  transmitted  signal.  Large  values  of  the 
correlation  integral  indicate  a  strong  resemblance,  while  small  values  indicate  weak 
resemblance.  The  values  of  the  time-delay  and  Doppler  shift  parameters,  which 
maximized  the  correlation  integral,  were  also  the  parameters  most  likely  to  represent  the 
target. 

Using  the  same  notion,  imaging  can  be  achieved  by  determining  the  position  and 
velocity  parameters  that  represent  the  best  “fit”  to  the  measured  data.  These  will  also  be 
the  parameters  that  localize  the  position  and  velocity  of  the  unknown  scatterers  in  phase 
space  (position-velocity)  space.  Cross  correlation  in  the  time  domain,  is  multiplication  in 
the  Fourier  transform  domain  (. Fourier  correlation  theorem)  -  hence  the  association  with 
filtering  since  filters  are  multiplication  operators  in  the  Fourier  transform  domain. 

The  model  developed  for  y/ mM(y,z,t)  depends  on  the  position  of  the  scatterer 
x  =  (xvx2,x3)  and  its  velocity  v  =  (v15v2,v3)  .  Using  this  form  for  y/scatt(y,z,t)  to  build  a 
six-parameter  scattering  model  based  on  the  arbitrary  position  p  =  ( prp2,p3)  and 
arbitrary  velocity  u  =  (ul,u2,u3),  the  objective  will  be  to  detennine  the  values  of p  and  u 
that  maximize  the  cross  correlation  between  the  arbitrary  model  and  the  measured  data. 

In  terms  of  the  Doppler  shift,  the  scattered  field  (from  known  scatterers  pv(x) 
with  positions  x  and  velocities  v)  is 

9  -ia>,(t-\z\/c) 

-CO  e  rr  C 

V scatty ’Z’0  *  ,  /  \2  i  n  ,  II  eXP  Tikv(P+  “  V({  ~  ZS(Z  ~  X)  1  C)]  } 

y\n )  |  z  ||  y  |  v  } 

x  ,s%(/  -  (|  z  |  -(p+  fi)gr)  /  c^pv(x)d'vd'x 
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The  effect  of  data  shifted  in  time  is  an  image  that  is  translated  in  range.  Since  this 
has  no  bearing  on  the  imaging  process,  the  scattered  field  result  can  be  simplified  by 
substituting  t'  =  t-\z\/c  which  yields 


Y(y,z,V). 


i  -icot' 

-co  e  y 
_ y _ 

(4 n)2  \z\\y 


■  JJ  exp  | -iky{y>+  ij)Sx  -  v(t'+  igc  /  c)\  J 

x  §/p(t '+  (j0+  ii)gr  /  c}pv(x)d3vd3x 


Ignoring  the  intensity  prefactor,  a  parametric  model  representing  the  field  from  an 
unknown  point  source  is 


y/(y,z,t')  =  -eUOyt  exp  &)§p  -  u(t'+  sqp  /  c)]  j 

x  §/o(t '+  (fi+  S)op  /  c  j 


(3.22) 


And  the  image  is  created  as  follows:  For  each  p  and  u,  determine  the  function 

I(p,u)  =  JJJ  y/scatt(y,z,t')y/Xy,z,t')dt'dmyd"z  (3.23) 

where  m,  n  depends  on  the  configuration  of  the  transmitter(s)/receiver(s),  and  the 
integrals  are  over  all  values  of  t\y  and  z  for  which  measured  data  is  available. 

I(p,u)  =  JJJ  i//scatt(y,z,t')i//\y,z,t')Q((o,t',p,u,y,z)dt'dmydnz 

I(p,u)  =  -JJJ 

X  §/o(t'+  (y>+  fi)g p  /  c)y/scatt (y, z,  t ')dt '  dmyd"z 
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IV.  IMAGING  POINT  SPREAD  FUNCTION 


This  chapter  begins  by  looking  at  the  radar  ambiguity  function,  and  the 
relationship  between  the  transmitted  waveform,  and  the  accuracy  with  which  target  range 
and  velocity  can  be  estimated.  The  imaging  point-spread  function  (PSF)  is  also  derived 
based  on  the  imaging  scheme  developed  in  Chapter  III. 

A.  RADAR  AMBIGUITY  FUNCTION 


For  a  stationary  target,  the  output  of  the  correlation  receiver  is  the  cross 
correlation  between  (1)  the  received  signal  plus  noise  and  (2)  the  transmitted  signal.  In 
many  radar  applications,  however,  the  target  is  moving  so  that  its  echo  signal  has  a 
Doppler  frequency  shift.  The  output  is  hence  a  cross  correlation  between  the  received 
signal  and  the  Doppler-shifted  transmitted  signal. 

The  nature  of  the  correlation  receiver  output  as  a  function  of  both  time  and 
Doppler  frequency  is  important  for  understanding  the  properties  of  a  radar  wavefonn,  in 
particular  its  effect  on  measurement  accuracy,  target  resolution,  and  ambiguities  in  range 
and  radial  velocity.  These  aspects  of  the  correlation  receiver  output  will  be  examined  in 
turn. 


When  the  received  echo  signal  is  large  compared  to  noise,  the  output  of  the 
correlation  receiver  may  be  written  as  the  following  (see  Equations  (2.15)  and  (2.16)) 


w  1 

ni-(v+v')(T-z') 

p(o',r')%(o  -  v',r  -  z')e  2  dr' do' 


where  %(o  -  o',  r  -  r')  is  the  radar  ambiguity  function.  The  problem  of  radar  imaging  is 
concerned  with  estimating  the  location  and  strength  of  the  point  scatterers  that  are 
assumed  to  make  up  the  target.  Consequently,  understanding  %  is  essential  as  the  end 
state  involves  ‘inverting’  the  radar  data  model.  By  definition,  we  have 

»  (  i  't  (  i  A 

Z(o-o',t-t')=  J  s.nc  ')js*nc  yt'+-{T-T')jeiiv~vy'dt'  (4.1) 
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1. 


Basic  Properties 


For  simplicity,  the  origin  is  chosen  to  be  the  true  target  time-delay  and  frequency- 
shift.  This  assumption  will  be  carried  throughout  the  rest  of  the  discussion.  Some  key 
properties  that  follow  from  Equation  (4.1)  [Refs.  14,  18  and  19]  are  listed 

Signal  energy  :  | %(u,r)\  <  |  j(0,0)|  (4.2) 

oo 

Ambiguity  volume  :  JJ  |^(n>,r)|  cl  id v  =  1  (normalized  signal)  (4.3) 

—oo 

Symmetry :  \%(-v,-t)\  =  \%(v,  r)\  (4.4) 


Frequency  -  domain  :  %(u,t) 


-uK  (co  +  -u)eimdco 

2  "!C  2 


(4.5) 


Equation  (4.2)  states  that  the  maximum  value  of  the  ambiguity  function  occurs  at  the 
origin  (true  target  time-delay  and  frequency-shift).  Equation  (4.3)  asserts  that  the  total 
volume  under  the  ambiguity  surface  is  a  constant.  This  gives  rise  to  a  radar  uncertainty 
principle :  choosing  a  signal  s(t)  so  that  the  ambiguity  surface  will  be  narrow  in  one 
dimension  will  cause  it  to  be  correspondingly  wide  in  the  other  dimension.  In  the  case  of 
range-Doppler  imaging,  a  waveform  with  good  range  resolution  has  poor  Doppler 
resolution  and  vice  versa. 


2.  Range  Resolution 


Equation  (4.5)  presents  an  intuitive  look  at  range  resolution.  Consider  a  fixed 
target  (v  =  0)  and  rewrite  Equation  (4.5)  as 


z(  o,0 


1 

2  n 


J  \s.nc((o)\2e-imd(o 


(4.6) 


Recognize  the  form  of  Equation  (4.6)  as  the  inverse  Fourier  transform  (of  the  transmitted 
power  spectral  density).  Then,  in  order  for  |^(0,r)|  to  have  high  resolution  in  z 
measurements  (ideally  a  delta  function),  the  power  spectral  density  has  to  be  identically 
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one;  such  a  signal  will  have  infinite  energy  or  broadly  supported  in  the  Fourier  domain. 
Hence,  it  is  clear  that  better  range  resolution  is  obtained  when  the  signal  bandwidth  is 
broad. 

A  more  rigorous  study  for  estimation  errors  (in  the  presence  of  Gaussian  noise)  is 
developed  in  [1]  by  examining  the  behavior  of  the  ambiguity  function  in  the 
neighborhood  of  its  main  peak.  This  study  shows  that  estimation  error  of  o  and  z  is 
described  by  an  ellipse,  and  that  time-domain  resolution  is  inversely  proportional  to 
frequency-domain  bandwidth.  Hence,  bandwidth  is  inherently  important  to  problems  in 
radar  target  imaging.  A  further  interesting  result  makes  use  of  the  definition  of 
bandwidth  in  terms  of  the  second  moments  (probabilistic  variance)  of  the  power  spectrum 
to  yield  (for  a  narrow-band  signal  s(t )  =  a(t)el<t,(n  ).  We  can  write 


da(t') 

dt' 


dt '  4 — 
2 


1  r( 

2  A 


(  dQ>(t')\ 


dt ’ 


a(t')  |2  dt '  -  Co] 


(4.7) 


where  fi  represents  bandwidth.  The  useful  relationship  illustrates  a  nonlinear  O(t) 
(phase  modulation)  will  increase  the  bandwidth  in  comparison  with  signals  that  are 
modulated  in  amplitude  only.  This  is  the  basis  for  pulse  compression  (see  Chapter  V),  a 
technique  of  practical  importance  since  it  allows  for  the  creation  of  fine  range-resolution 
waveforms  that  are  also  of  long  duration  (with  sufficient  energy  on  target). 


3.  Doppler  Resolution 


To  study  Doppler  resolution,  consider  a  target  whose  range  is  known.  Thus, 
setting  r  =  0  in  Equation  (4.1)  yields 


00 

\z(»>Q)  |= 

J  \2e‘u,'dt' 

-00 

(4.8) 


It  is  evident  that  better  Doppler  resolution  is  obtained  from  a  long  duration  signal.  This  is 
in  stark  contrast  to  range  resolution  requirements.  In  addition,  the  Doppler  resolution  is 
determined  only  by  the  amplitude  modulation  of  the  signal  and  not  by  the  phase. 
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B. 


IMAGE  ANALYSIS 


From  Chapter  III,  an  explicit  form  for  the  image  is 

J I  \  fff  ,,  S  iaj'  ik Ap+S)6p-u(t'+iQplc)\ 

I(p,u)  =  -\\\Q(a>,t',p,u,y,z)e  1  e  y 

(4-9) 

x  &°y  '+  (J0+  ®)<3P  /  c)i//scatl(y,z,t')dt'dmyd"z 

To  examine  the  performance  of  the  imaging  scheme,  the  ideal  (expected)  data  are 
inserted  back  into  the  image  equation.  Substituting  for  y/ scan(y,z,t')  yields 


T,  .  rrr  rrC02Q(CO,t',p,U,y,z) 

=  JJJ  JJ  '  -  .1 1  I,  | — pSx) 

x  exp  ^k  (y>+3)Qp  -x  —  (u  —  v)t'—  u(£Qp)  /  c  +  v  (rg t)/c]} 
x  .%([ '+  (j0+  S)gp  /  c '+  (j0+  i))gc  /  c yi 3 vd \xdl ' d yd " z 


(4.10) 


The  distances  |  z  |  and  |  y  \  are  scaling  amplitude  factors  and  not  essential  for 
imaging  purposes.  Thus,  choose 

(4;r)2  I  z  II  v  I 

Q(co,t\p,u,y,z)  = - - — — J(p,u,y,z )  (4.11) 

a> 

y 


where  J(p,u,y,z)  depends  on  the  geometry  and  is  chosen  to  compensate  for  a  Jacobian 
that  results  from  the  integral  of  the  variable  t  ’ 

I(p,u )  =  JJJ  JJ.?/o(t'+  (j0+  i))gn  /  c)%(f+  (;d+  i))gr  /  c)e  ,ky{y+z)Au-v)t 

x  exp  ^ky(p+  i5)c[/?  - x  -  u(xap)  /  c  +  v(rgr)  /  c]  J-  (4. 12) 

x  J(p,u,y,z)pv(x)dt' d3vd3xdmydnz 

which  is 

I(p,u)  =  JJ  K(p,u,y  ,z)pv(x)d\d3  x  (4.13) 

where 
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K(p,u,y,z )  =  JJ  exp  { -iky(y>+  £))$u{£qp)  -  v(igr)]  /  c  j 
x  | ,%(/ '+  (j0+  n)gp  /  c)%(l  '+  (y>+  T>)gr  /  !^<s 

x  J(p,u,y,z)d"'ydnz 


-ik  (p+z>)^u-v)t' 


dt' 


(4.14) 


is  the  point-spread  function  describing  the  behavior  of  the  imaging  system.  Apply  the 
change  of  variables  t  =  t'+  ~(P+  ®)&P  +  x)  /  c  and  set  r  =  (y>+  i))c(  p  -  x)  /  c  to  obtain 


K(p,u,y,z)  =  JJ  exp  { -iky(y>+  B)$u(xgp)  -  v(igr)]  /  c  j 

x  exp  jz'k;(j0+  i>)<J(p-  x)  +  ^(m  -  v)(o+  £fap  +  x)]  /  cj 


(4.15) 


f  nJ  1  ^  n/f  1^  -ik  (ii>+z>)tfii-v)t  r 

J  %yt  +  -z j%yt - -z je  v  a  dt 


J(p,u,y,z)d"‘ydnz 


The  integral  in  square  brackets  is  the  radar  ambiguity  function 


»  f  1  A  f  1  'N 
z(v,  r)  =  J  ^  t  +  -  rj  .?/o^  t  -  -  Tj  e~iv,dt 


X  is  determined  by  the  transmitted  signal  and  |  %  \  attains  its  maximum  value  at  p  =  x , 
u  =  v  (for  general  configurations  of  transmitter/receiver).  So  the  imaging  PSF  can  be 
written  as 

K(p,u,y,z )  =  JJ  exp  ^-iky(j»+ fi)Qu(£gp)  -  v(ijgr)]  /  c  j 

x  exp  j/7cv(fi+  i>)<J(p-  x)  +  ^(u-  v)(p+Z))$p  +  x)]  /  c|  (4. 16) 

x  x(ky(?+  ®)9A  -  v),(;p+  £&p  - x) /  c^y (p ,u,y ,z)dm ydn z 
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V.  PULSE  COMPRESSION 


Detection  of  signals  in  noise  and  resolving  capability  are  both  fundamental  to 
radar  imaging.  As  discussed  in  Chapter  II  and  IV,  the  use  of  correlation  reception  and  a 
transmitted  waveform  with  non-linear  phase  modulation  provides  means  to  enhance 
performance  in  detection  and  resolution.  Pulse  compression  is  a  technique  that  fuses 
these  concepts  and  involves  transmitting  long  coded  pulses,  together  with  signal 
processing  methods  such  as  correlation  reception.  This  is  of  practical  importance  since  it 
allows  for  the  creation  of  fine  range-resolution  waveforms  that  are  also  of  long  duration 
(this  increases  radiation  energy  directed  on  target  without  relying  solely  on  boosting  radar 
transmitted  power).  This  chapter  gives  a  physical  interpretation  of  pulse  compression 
and  looks  at  a  specific  phase-coded  signal  known  as  the  chirp. 

A.  REVISITING  CORRELATION  RECEPTION 

Recall  in  Chapter  II,  the  output  of  the  correlation  receiver  is  a  correlation  integral 
between  the  received  signal  and  a  time-delayed,  Doppler-shifted  version  of  the 
transmitted  signal.  It  is  illustrative  to  explain  the  phenomenon  of  pulse  compression  with 
a  graphical  example.  Using  the  simple  one-dimensional  imaging  case  in  Equation  (2.12), 
the  correlation  receiver  output  can  be  written  as 

co 

7(0  =  J  0(Ov«,0'+  f)dt'  (5.4) 

-co 

which  is  a  correlation  between  5  and  5  .  If  5  =  5  (reflectivity  function  is  1),  Equation 

(5.4)  is  called  an  autocorrelation.  In  Figure  5.  the  signal  at  the  top  represents  the 
transmitted  waveform.  The  lower  graphs  represent  the  received  signal  for  different  shifts, 
according  to  Equation  (5.4).  When  a  signal  shifted  by  t  is  multiplied  by  the  waveform  at 
the  top  and  the  product  is  integrated  over  t  ’,  the  resulting  number  is  plotted  at  position  t 
on  the  graph  to  the  right. 
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s(t’) 


Figure  5.  Phenomenon  of  pulse  compression  where  the  energy  is  concentrated  at  a 

single  delay  time  [18]. 

Figure  5  shows  that,  although  the  output  signal  has  a  longer  time  duration  than  the 
original  signal  s,  its  energy  is  more  concentrated  at  a  single  delay  time.  It  is  a 
consequence  of  the  Cauchy- Schwarz  inequality  that  the  highest  peak  occurs  when  the 
signals  are  not  shifted  relative  to  each  other.  The  smaller  peaks,  called  sidelobes,  are 
undesirable  and  an  open  problem  is  to  design  pulses  whose  autocorrelations  have  the 
lowest  and  fewest  sidelobes.  An  issue  that  has  direct  bearing  on  signal  detection  as  the 
received  signal  is  expected  to  “look”  like  a  time-delayed,  Doppler-shifted  version  of 
itself. 


B. 


PHASE  CODING  -  CHIRPS 


A  radar  system  makes  the  most  efficient  use  of  power  when  it  transmits  a 
constant-amplitude  waveform.  Consequently,  most  radar  waveforms  are  designed  with 
variations  only  in  the  phase  of  the  signal.  The  most  important  and  most  commonly-used 
radar  wavefonn  is  the  chirp.  A  chirp  is  a  constant-amplitude  signal  whose  instantaneous 
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frequency  varies  linearly  with  time;  chirps  are  therefore  also  called  Linearly  Frequency 
Modulated  (LFM)  waveforms.  Linear  variation  of  the  instantaneous  frequency  implies 
that  d(j)  /  dt  =  comm  +  yt .  The  coefficient  y  is  called  the  (angular)  chirp  rate.  A  chirp  is 
thus  of  the  form 


s(t )  =  e 


i(co  .  t+—yt2 ) 


u[o  ,r]W 


where  w  (?)  is  1  in  the  interval  [0,  7]  and  0  otherwise. 

A  chirp  with  positive  chirp  slope  is  called  an  upchirp;  one  with  a  negative  chirp 
slope  is  a  downchirp.  The  instantaneous  (angular)  frequency  of  the  chirp  varies  from 
co  to  co  =  co  .  +  yt ,  with  a  center  frequency  of  <on  =  co  .  +  yt  /  2  .  Thus  the 

instantaneous  frequency  of  the  chirp  varies  over  the  interval  yco0-Ty  /  2, a>0  +  Ty  /  2  j. 
The  power  spectrum  of  a  chirp  is  roughly  constant  over  the  frequency  band 
[<z>0  -  Ty  /  2,co a  +  Ty  /  2] ,  thus  the  angular  frequency  bandwidth  Ty  can  be  estimated 
by  looking  at  the  range  of  instantaneous  frequencies  (see  Figure  6). 


iWj  -  r-f/1  .Ji.i  i  ry/2 


Figure  6.  Spectrum  for  a  chirp  pulse  of  length  r  [1]. 

Chirps  are  very  commonly  used  in  radar  imaging.  Pulse  compression  for  chirps 
has  the  following  intuitive  interpretation:  When  the  pulse  is  transmitted,  different  parts  of 
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the  pulse  are  coded  with  different  instantaneous  frequencies.  In  the  pulse  compression 
process,  the  different  frequencies  are  delayed  by  different  amounts  of  time,  so  that  all  the 
energy  emerges  at  the  same  time  [18]. 

C.  AMBIGUITY  FUNCTION  FOR  COMMON  WAVEFORMS 


From  the  preceding  discussion,  the  “ideal”  ambiguity  function  will  typically 
consist  of  a  single  peak  of  infinitesimal  thickness  at  the  origin  and  be  zero  everywhere 
else.  It  will  closely  approximate  an  impulse  function,  and  have  no  ambiguities  in  range 
or  Doppler  frequency.  The  infinitesimal  thickness  at  the  origin  will  pennit  the  time  delay 
and/or  frequency  to  be  determined  simultaneously  to  as  high  a  degree  of  accuracy  as 
desired.  Such  a  highly  desirable  ambiguity  function,  however,  is  not  theoretically 
allowed.  The  properties  of  the  ambiguity  function  implicitly  mandate  that  trade-offs  exist 
between  the  measured  dimensions. 

In  practice,  other  than  the  trade-offs  in  dimension  resolution,  waveforms  also 
generally  have  significant  response  outside  the  narrow  region  in  the  near  vicinity  of  the 
origin,  i.e.  there  are  many  possible  ambiguous  (blind)  range  and  Doppler.  Ambiguities 
are  a  consequence  of  discontinuous  waveforms,  such  as  in  a  pulse  train.  The  resolution 
for  single  rectangular  pulse,  single  chirp  and  coherent  pulse  trains  will  be  discussed  next. 


1.  Resolution  for  a  Single  Rectangular  Pulse 


Time  domain  impulses  have  infinite  bandwidth,  hence  they  are  idealized  signals. 
More  realistic  wavefonns  can  be  modeled,  for  example,  as 


Jl  if  0  <  t  <  T, 

[0  otherwise. 


(5.6) 


for  which  s  (t  -  1  /  2 t)s*  (t  +  1  /  2r)  is  non-zero  in  the  interval  t  e(|  r  |  /2,  T-  |  r  |  /2) . 
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The  ambiguity  function  is  then 


r  1  1  •  ,  rr-M/2 

z(u, t)=  [  j.  (*'--?>:  (t'+-z)ew‘ dt'  =  f 

/tV  ’  '  J  mcv  9  '  9  '  J|r[/2 


=  — e 
iu 


iv(T/2) 


eiv(TI2-\r\l2)  _  e-iu(TI2-\z\l2) 


(5.7) 


=  (r-  I  r  |)e'u772sinc 


1 


■u(T-|  r|) 


For  o  =  0 ,  then  the  delay  resolution  (by  the  peak-to-first-null  measure)  as 
observed  from  the  ambiguity  surface  |^(0,r)|2,  is  A r  =  T  (zero-crossing  at  t  =  T), 
which  corresponds  to  a  range  resolution  of  A R  =  cT  /  2.  For  r  =  0 ,  the  ambiguity 
surface  |  %(v, 0)  |2  yields  Au  =  1/7’  since  the  first  zero  of  the  sine  function  occurs  when 
its  argument  is  n .  The  Doppler  shift  is  related  to  down-range  relative  velocity  by 
vD  =  -2 vuQ  /  c  (see  Equation  (2.4)),  and  the  corresponding  velocity  resolution  is 
Av  =  A0  /  IT . 


2,  Resolution  for  a  Single  Chirp 


This  is  another  example  of  a  ridge  type  ambiguity  function  produced  by  linearly 
frequency  modulating  a  rectangular  pulse  over  a  bandwidth  j3 .  The  pulse  width  T  is 
large  compared  to  1  /  /? .  The  frequency  modulation  increases  the  spectral  bandwidth  of 
the  pulse  so  that  J3T  ?  1 .  Since  the  pulse  width  T  and  the  bandwidth  [J  can  be  chosen 
independent  of  one  another,  the  time-delay  and  frequency  accuracies  are  independent  of 
the  other. 


The  chirp  signal  is  one  whose  frequency  changes  linearly  with  time 


r  i  ^ 


0 


if  0  <  t  <  T, 
otherwise. 


(5.8) 
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The  constant  y  is  known  as  the  chirp  rate.  The  ambiguity  function  is  then 


OO-I  < 

r  1  I  rT-W  1/2 

%(v,z)  =  [  j.  (t'--z)st  (t'+-T)ew‘ dt'=  f  e,iv+n)t 

v  7  7  J  incK  y  /  mcy  9  7  Jiri/2 


dt' 


1 


i(u+rr)TI2 


i(v  +  yz) 

=  (T-\t  \)e,{v+rT)T,2smc 


ei(v+rrXT/2-\r\/2)  _  /2-\v\/2) 


(5.9) 


^(u  +  yr)(T-\T\) 


Observe  that  a  chirp  leads  to  an  ambiguity  function  that  is  a  smc-function  in  both  the  time 
and  frequency  domain.  It  is  instructive  to  examine  this  property  again  by  looking  at  the 
signal  resolution. 

For  r  =  0 ,  similar  result  of  At>  =  1  /  T  is  obtained  as  in  the  case  of  the  rectangular 
pulse.  This  is  because  the  amplitude  modulation  is  identical  in  both  cases. 
Consequently,  Doppler  resolution  is  the  same.  However,  with  phase  modulation,  it  is 
expected  that  the  delay  resolution  (hence  range  resolution)  will  be  enhanced. 

For  v  =  0 ,  when  r  >  0  the  first  null  occurs  when  n  =  {T  -  r)yr  /  2 .  Reference 
[19]  shows  that  the  expression  n-(T  -  r)yv  /  2  has  a  quadratic  equation  form  with  a 
solution  /S.z  =  2n I  yT ,  which  corresponds  to  a  range  resolution  of  AR  =  c/  fd  where 
/?  =  2 nyT  is  the  bandwidth  (in  Hertz).  For  r  =  0,  Ao  =|  /(^,0)  |2=  1  /  T  since  the  first 
zero  of  the  sine  function  occurs  when  its  argument  is  n .  The  Doppler  shift  is  related  to 
down-range  relative  velocity  by  uf)  =  -2vt;0  /  c  (see  Equation  (2.4)),  and  the 

corresponding  velocity  resolution  Ay  =  A()  /  2T . 

The  pulse  compression  ratio  is  defined  as  the  product  of  pulse  spectral  bandwidth 
/ 3  and  the  uncompressed  pulsewidth  T.  Comparing  the  range  resolution  between  the 
single  rectangular  pulse  and  single  chirp, 


A R 


rectangular 


A R 


chirp 


cT  /  2  pT 

dp  ~  2 


(5.10) 
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In  other  words,  phase  modulation  improves  the  range  resolution  by  a  factor  of  one-half 
the  time-bandwidth  product. 

3.  Resolution  for  Coherent  Pulse  Trains 


Following  the  analysis  given  in  [18],  a  train  of  N  identical  pulses  is  considered: 


1  JV— 1 

S(t)  =  ~  nTR) 

yj  N  n= 0 


(5.11) 


where  T  is  the  pulse  repetition  interval,  which  might  be  on  the  order  of  1ms.  The 
reciprocal  1  /  TR  is  called  the  pulse  repetition  frequency  (PRF). 


The  ambiguity  function  for  the  pulse  train  is 


1  iV-lJV-l 

n=0  m=0 


u*(t  +  t'-  mT  nT  )eWT  dr' 


1 

N 


N- 1  N-\ 

X  X  xu  (t  -  (n  -  m)TR » 

n= 0  m= 0 


(5.12) 


where  xu  denotes  the  ambiguity  function  for  it  and  where  we  have  used  the  substitution 
t"  =  r-  nT  .  A  lengthy  calculation  involving  rearrangement  of  the  order  of  summation 
and  the  summing  of  a  geometric  series  results  in 


1  JV— 1 

X(o,t)  =  —  X  Xu(r-pTR,o)e 

^  p=-(N-\) 


iuTR(N-\-\p\) 


sin  ttuTr(N-  |  p  |) 
s\n(7wTR ) 


(5.13) 


If  the  pulses  u  are  sufficiently  well-separated  so  their  ambiguity  functions  %  do  not 
overlap,  then  we  take  the  absolute  value  of  Equation  (5.13)  to  obtain 


T) 


i  N-\ 

~  X 

NPJT»- 1) 


Xu(t-PTr,v) 


sin  ttuTr{N-  |  p  |) 
sin  (ttuTr) 


(5.14) 
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Such  an  ambiguity  function  has  a  “bed  of  nails”  appearance  (see  Figure  7),  with  peaks  at 
r  =  pTR,  />  =  -(Af-l),-(7V-2),...,0,l,...,(7V-l)  and  o  such  that 
ttvTrN  =  (m  +  1  /  1)n ,  m  =  0,1,2,... 


Figure  7.  Ambiguity  function  for  a  train  of  pulses. 


Following  previous  analysis,  consider  the  u  =  0  cut  of  the  ambiguity  surface  to 
estimate  the  range  resolution  yields 

l*(0,r)H-j-  Z  \xM-PTRm(N-\p\)  (5.15) 

^  p=-(N-l) 


where  we  have  used  the  fact  that  lim[sin^x]  /  sin  x-b.  We  see  that  the  range  resolution 

x-»0 

of  the  main  peak  is  the  same  as  the  range  resolution  of  the  individual  pulse  u,  but  that  we 
now  have  range  ambiguities  due  to  the  extraneous  peaks  at  t  =  pTR  . 


To  estimate  the  Doppler  resolution,  consider  the  r  =  0  cut  of  Equation  (5.14)  and 
for  well-separated  pulses,  only  the  p=0  term  contributes 


i  JV-l 
™  p=-{N- 1) 


sinf^W] 
sin  [uTr] 


(5.15) 
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Consider  the  numerator  argument  since  it  varies  much  faster.  It  shows  that  the  Doppler 
resolution  is  1  /  NT  ,  and  hence,  the  Doppler  resolution  of  the  pulse  train  is  higher  than 
that  of  a  single  pulse. 

The  many  ambiguities  produced  by  a  pulse  train  may  appear  to  lead  to  a  poor 
radar  wavefonn.  The  time-delay  measurement  accuracy  (detennined  by  pulsewidth  T) 
and  the  frequency  accuracy  (detennined  by  the  pulse  repetition  interval  TR)  can  be 

selected  independently.  If  pulse  repetition  interval  T  is  such  that  no  radar  echoes  are 
expected  with  a  time  delay  greater  than  T  ,  and  no  Doppler-frequency  shifts  are  expected 
greater  than  1  IT,  then  the  effective  ambiguity  function  reduces  to  just  a  single  spike  at 
the  origin  whose  dimensions  are  detennined  by  T  and  TR .  In  practice,  ambiguities  can 
also  be  resolved  with  different  PRF,  and  many  radars  employ  this  type  of  waveform. 
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VI.  RESULTS  AND  CONCLUSION 


Here,  analysis  of  the  imaging  point-spread  function  (PSF)  (Equation  (4.16))  is 
carried  out  using  Matlab.  Equation  (4.16)  does  not  restrict  the  choice  of  the  transmitted 
wavefonn  nor  the  transmitter/receiver  configuration.  The  goal  here  is  to  study  the 
behavior  of  the  PSF  under  different  waveforms  and  transmitter/receiver  configuration. 

Figure  8  shows  the  transmitter/receiver  geometries  for  the  simulation  performed. 
The  receivers  are  denoted  by  “zl”  to  “z8”  while  the  transmitter  is  denoted  by  “yl”.  For 
configurations  of  two  to  four  receivers,  the  geometry  as  shown  in  Figure  8(a)  is  used, 
while  Figure  8(b)  shows  the  geometry  for  eight  receivers. 


x-component 

(position) 


(a) 
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x-component 

(position) 


(b) 

Figure  8.  Transmitter/Receiver  configuration  for  2,  3,  4  and  8  receivers. 

A.  IMAGING  PSF  FOR  SINGLE  RECTANGULAR  PULSE 

The  transmitted  waveform  is  single  rectangular  pulse  defined  in  Equation  (5.6). 
Its  associated  ambiguity  function  (Equation  (5.7))  is  substituted  into  the  imaging  PSF, 
Equation  (4.16)  and  implemented  in  Matlab.  Equation  (4.16)  is  capable  of  determining  a 
three-dimensional  distribution  of  positions  in  space  and  the  corresponding  three- 
dimensional  velocities  (i.e.  a  six-dimensional  phase  space  image).  However,  to  study  the 
effects  of  transmitter/receiver  geometry  and  translation  in  position/velocity  on  the 
imaging  PSF,  separate  two-dimensional  plots  of  the  PSF  in  position  and  velocity  space 
will  suffice. 
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1. 


Localization  of  Point-Scatterer  in  Position  Space 


The  point  scatterer  is  located  at  position  (xt,x2)  =  (0,0).  Two-dimensional  plots 

of  the  imaging  PSF  in  position  space  for  two,  three,  four  and  eight  receivers  are  given  in 
Figures  9(a)-(d)  respectively.  Clearly,  the  point  scatterers  are  localized  in  position  space 
about  the  scatterer  position  (0,0).  However,  under  the  different  geometries,  the  point 
scatterer  is  represented  differently,  with  the  ‘two  receivers’  case  resembling  the  least  like 
a  point.  In  Figure  9(a),  the  point  scatterer  appears  longish  and  rectangular  in  shape. 
These  distortions  of  the  point  scatterer  are  responsible  for  the  image  artifacts  encountered 
during  imaging.  The  extent  of  the  PSF  distortions  will  detennine  the  level  of  artifacts 
produced.  In  Figure  9(d),  the  PSF  appears  as  a  blob  and  one  will  expect  to  see  smudges 
in  the  resulting  image.  It  can  be  observed  that  the  geometry  of  the  transmitters/receiver 
plays  a  role  in  shaping  the  PSF. 


(a) 


(b) 
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x-component  x-component 

(velocity)  (velocity) 


(c)  (d) 

Figure  9.  2-D  Localization  of  point-scatterer  in  position  space  (for  rectangular  pulse). 

Figure  10  shows  the  point-scatterer  is  relocated  to  position  (xvx2)  =  (5,5)  .  The 

resulting  imaging  PSF  becomes  localized  now  at  the  new  position  (5,5).  The  PSF  retains 
its  shape,  and  the  characteristic  of  translation  (in  position)  invariance  for  the  PSF  is 
necessary  for  a  good  imaging  system. 


x-component 

(velocity) 


Figure  10.  Position  invariance  of  imaging  PSF. 

2.  Localization  of  Point-Scatterer  in  Velocity  Space 

The  point  scatterer  is  located  at  velocity  (v1?v2)  =  (0,0) .  Two-dimensional  plots 
of  the  imaging  PSF  in  velocity  space  for  two,  three,  four  and  eight  receivers  are  given  in 
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Figures  ll(a)-(d)  respectively.  Clearly,  the  point  scatterers  are  also  localized  in  velocity 
space,  and  similar  to  the  case  in  position  space,  the  geometry  of  the  transmitters/receiver 
plays  a  role  in  the  shape  of  the  PSF.  The  velocity  resolution  is  visibly  poorer. 


Figure  1 1 .  2-D  Localization  of  point-scatterer  in  velocity  space  (for  rectangular  pulse). 

When  the  point-scatterer  is  moving  with  a  new  velocity  of  (v1?v2)  =  (-35,-50) , 

Figure  12  shows  that  the  resulting  imaging  PSF  becomes  localized  now  at  the  new 
velocity  (-35,-50) .  The  imaging  PSF  demonstrates  similar  translation  invariance  in  the 
velocity  space. 
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x-component 

(velocity) 


Figure  12.  Velocity  invariance  of  imaging  PSF. 

B.  IMAGING  PSF  FOR  SINGLE  CHIRP 

In  Chapter  V,  it  was  shown  that  a  chirp  signal  improves  range  resolution.  The 
associated  ambiguity  function  for  a  single  chirp  (Equation  (5.9))  is  substituted  into  the 
imaging  PSF,  Equation  (4.16)  and  implemented  in  Matlab. 

1.  Localization  of  Point-Scatterer  in  Position  Space 

Again,  the  point  scatterer  is  located  at  position  (x^x,)  =  (0,0) .  Two-dimensional 
plots  of  the  imaging  PSF  in  position  space  for  two,  three,  four  and  eight  receivers  are 
given  in  Figures  13(a)-(d)  respectively.  The  PSF  is  visibly  better  localized  as  compared 
to  transmitting  a  rectangular  pulse. 
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Figure  13.  2-D  Localization  of  point-scatterer  in  position  space  (for  chirp  signal). 


C.  CONCLUSION 

Image  artifacts  result  when  motion  in  the  imaging  scene  is  not  addressed.  In 
recent  years,  a  number  of  attempts  to  develop  imaging  techniques  that  can  handle  moving 
objects  have  been  proposed,  however,  all  of  these  techniques  rely  on  making  the 
approximation  that  a  target  in  motion  is  assumed  to  be  momentarily  stationary  while  it  is 
being  interrogated  by  a  radar  pulse.  This  thesis  is  a  study  of  the  new  imaging  approach 
developed  by  Cheney,  Borden  [2]  that  can  accommodate  target  motion  during  the 
imaging  process. 

In  this  thesis,  the  subject  of  radar  imaging  from  scattered  waves  is  explored  and 
applied  to  derive  the  corresponding  imaging  point-spread  function  for  the  new  imaging 
approach.  The  simulation  result  (using  Matlab)  obtained  showed  that  the  imaging  point- 
spread  function  is  well  behaved,  specifically  it  localizes  the  target  in  phase  space  and  is 
translation  invariant.  In  addition,  it  is  apparent  that  the  geometry  of  the  transmitters  and 
receivers  disposition  affects  the  behavior  of  the  point-spread  function,  and  hence  the 
imaging  system. 

This  work  represents  the  early  stages  of  development  of  a  full  imaging  algorithm 
applicable  to  moving  targets  in  a  multistatic  environment.  Follow-on  work  will  include 

further  simulation  efforts.  This  thesis  used  single  pulse  radar  waveforms  and  a  single 
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transmitter,  but  in  reality,  fielded  radar  systems  typically  utilize  pulse  compression 
techniques  and  send  out  long,  coded  pulse  trains.  Neither  need  the  multistatic  radar 
system  be  restricted  to  one  transmitter.  In  this  case,  the  imaging  point-spread  function 
has  to  be  adjusted  to  incorporate  more  realistic  radar  waveforms  and  optimize  the 
transmitter/receiver  geometries  for  implementing  the  imaging  algorithm.  These  will  be 
used  to  support  the  development  of  the  eventual  imaging  algorithm.  Finally,  real  world 
target  data  can  then  be  applied  to  the  developed  imaging  algorithm  to  assess  its 
performance. 
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APPENDIX:  MATLAB  CODES 


%%%%%%%%%%%%%%%%% 

%  This  program  evaluates  the  Point  Spread  Function  (PSF) 

%  The  program  uses  the  approach,  techniques  on  imaging  from  Professor 
Borden,  %  Brett 

%  It  is  a  six-para  me  ter  scattering  model  based  on  arbitrary  position 
%  p=(px,py,pz)  &  velocity  u=(ux,uy,uz) 

%%%%%%%%%%%%%%%%% 


clear  all 


%%  Variables  to  be  s  p  e  c  i  f  i 

x  =(  10A3) *[ 0;  0;  0] ; 
v  =[  0;  0;  0] ; 
y 1 =[  2;  0;  0] ; 
z  1  =[  -  5;  5;  0]  ; 
z 2  =[  15; 15;  0] ; 
z 3  =[ -  2;  -  15;  0]  ; 
z 4  =[ -  10;  10;  0]  ; 
z 5  =[  0;  14.  1;  0]  ; 
z 6  =[  10;  10;  0]  ; 
z 7  =[  14.  1;  0;  0]  ; 
z  8  =[  10;  -  10;  0]  ; 
u=[  0;  0;  0]  ; 
pul  se-10A(  -  4)  ; 


%  a  c  t  u  a  I  target 
%  a  c  t  u  a  I  target  velocity 
%  t  r  a  n  s  mi  1 1  e  r  _  1 
%  r  e  c  e  i  v  e  r  _  1 
%  r  e  c  e  i  v  e  r  _  2 
%  r  e  c  e  i  v  e  r  3 


%%  P  a  r  a  me  t  e  r  s  Defined 

Y 1  =  (  1/  sqrt  (  y  1  (  1)  A  2  +  y  1  (  2)  A  2  +  y  1  (  3)  A2)  )  .  *  y  1 ; 

Zl=(  l/sqrt(zl(  1)  A2+zl(  2)  A2+zl(  3)  A2)  )  .  *zl; 

Z  2  =  (  1/  sqrt  (  z  2  (  1)  A  2  +  z  2  (  2)  A  2  +  z  2  (  3)  A2)  )  .  *  z  2 ; 

Z  3  =(  1/  sqrt  (  z  3  (  1) A2  +z  3 ( 2) A  2  +z  3 ( 3) A2) ) .  *  z  3 ; 

Z4  =  (  1/  sqr  t ( z  4 ( 1) A  2  +  z  4 ( 2) A  2  +  z  4 ( 3)  A2)  )  .  *  z  4 ; 

Z5  =  (  1/  sqrt  (  z  5  (  1) A  2  +z  5 ( 2) A2  +z  5 ( 3) A2) ) . *  z  5 ; 

Z6  =  (  1/  sqrt  (  z  6  (  1)  A  2  +  z  6  (  2)  A  2  +  z  6  (  3)  A2)  )  .  *  z  6 ; 

Z  7  =  (  1/  sqrt  (  z  7  (  1) A  2  +  z  7  (  2) A  2  +  z  7 ( 3) A2) ) ,  *  z  7 ; 

Z  8  =  (  1/  sqrt  (  z  8  (  1) A  2  +  z  8  (  2) A  2  +  z  8 ( 3) A2) ) ,  *  z  8 ; 

J  =  1 ;  %  J  a  c  o  b  i  a  n 

c=3* 10 A8; 

k_y-2*pi  *( 10A10) / c;  %  wave  vector  @f  = 1 0  GHz 


%%  PSF  Range  Defined 

p  x  st  art  =- 20; 
p_  x_  e  n  d  =2  0 ; 
p'y'st  art  =- 20; 
p_y_end=20; 

res_p=0.  1;  %  r es o I  u t i  o n  o f  p 

p_x=p_x_start :  res_p:  p_x_end;  %  define  vector  of  x-axis  values 

P_y-P_y_st art :  res_p;  p_y_end;  %  define  vector  of  y-axis  values 

%%  PSF  Velocity  Defined 

u  x  start  —  2  0  0; 
u_x_end=200; 
u  _  y  _  s  t art  =- 150; 
u_y_end=150; 
r  e  s  _  u  =1 ; 

u_x=u_x_start :  r  e  s  _  u :  u  _  x  _  e  n  d ; 
u.y=u_y_st  art :  r  e  s  _  u :  u  _  y  _  e  n  d ; 
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Localization  in  Position-Space 


nx=length(p_x);  %  length  of  x  =  columns  in  K 

ny=length(p_y);  %  length  of  y  =  rows  in  K 

K=zeros(ny,nx);  %  initialize  K  matrix 

Kl=zeros( ny,  nx)  ; 

K2=zeros( ny,  nx)  ; 

K3=zeros( nx,  ny)  ; 

K4  =z  e  r  o  s  (  n  x ,  n  y )  ; 

K5  =z  e  r  o  s  ( nx,  ny)  ; 

K6  =z  e  r  o  s  ( nx,  ny)  ; 

K7=zeros( nx,  ny)  ; 

K8=zeros( nx,  ny)  ; 

for  m=l :  n  y 

for  n  =1 :  n  x 

p=(10A3)*[p_x_start+((n-l)*res_p);p_y_end-((m-l)*res_p);0] 

%  f  o  r  receiver  1 

a  r  g  1  1=(  (  Yl+Zlj.  1  ) *  ( (  u  * ( ( Z 1 . 1 )*p))-(v*((Zl. 1 ) *  x ) ) ) ; 
phasel_l=exp( - i  *k  y*argl_l/ c) ; 

a  r  g  2  1=(  (  Yl+Zl)  .  1  T*((p-xy+(0.  5*  (  u  -  v)  *  (  (  (  Y1 +Z1 )  .  1  ) *( P+x) ) ) ) 
p  h  a  s  e  2  1  =e  x  p  (  i  * k  y  *  a  r  g  2  1/ c)  ; 

t  a  u_  1  =1/  c  M  (  (  Y 1  +Z1 )  .  1  )  *  [  p-  x )  )  ; 

ne u_ 1 =k_  y *  (  (  (  Y 1  +Z1 )  ,  1  ) * (  u -  v )  )  ; 

%r  e  c  t  a  n  g  u  I  a  r  pulse 

a  mb  _ 1 =p  i  *  ( pul  se-abs(tau  1 ) ) *  e  x  p ( - 
i  *0.  5*neu_l*pui  s  e ) *  s i  n  c ( ( 0 .  5/pi  J*( - neul) *(  pul  se-abs(tau_l))); 

Kl(  m,  n)  =phasel_l*phase2_l*amb_l*J  ; 

%  f  o  r  receiver  2 

a  r  g  1  2 =(  (  Y 1  +Z2  T .  1  )  *  (  (  u *(  (  Z2 .  1  )* p )  )  -  (  v *  (  (  Z2 .  1  )  * x )  )  )  ; 
ph  a  s  e 1_  2  =ex  p( - i  *  k  y  *  a  r  g  1  2/c); 

a  r  g  2  2  =(  (Yl+Zl)  .  1  I*(  (  p  -  xJ+(  0.  5*(u-v)*(((Yl+Z2).  1  ) *( p+x) ) ) ) 
p  h  a  s  e  2  2=exp(  i  *  k  y  *  a  r  g  2  2/  c)  ; 
t  a  u_  2  =1/  c  M  (  (  Y 1  +Z2 )  .  1  )  *  [  p-  x )  )  ; 
neu"2=k  y*( ( ( Y1+Z2) . 1 ) *( u- v)  )  ; 
amb  2=pF*(  pul  se  -  abs(tau  2 ) )  *  e  x  p  ( - 
i  *0.  5*neu_2?pul  s  e ) *  s i  n  c ( ( 0 .  5/pi J*( - neu_2) *( pul  se-abs(tau_2))); 

K2(  m,  n)  =phasel_2*phase2_2*amb_2*j  ; 

%  f  o  r  receiver  3 

a  r  g  1  3=(  (  Y1+Z3T.  1  )  *  ( (  u  *  (  (  Z  3 .  1  )*p))-(v*((Z3.  '  )  *  x )  )  )  ; 
phasel  3  =e  x  p ( - i *  k  y  *  a  r  g  1  3/c); 

a  r  g  2  3=(  (  Y1+Z3)  .  ,T*((p-xT+(0.  5*(  u-  v)  *(  (  (Y1+Z3)  .  '  )  *(  p+x)  )  )  ) 
phase2_3=exp(  i  *k_y5|:arg2  3/c); 
t  au  3=1/  cM  (  (  Yl+23)  .  1  )  *[p-  x)  )  ; 
neu~3=k  y * ( ( ( Y1+Z3) .  1 ) *( u- v)  )  ; 
amb_3=pf*(  pul  se-abs(tau  3 ) )  *  e  x  p  ( - 
i  *  0 .  5  *  n  e  u  3* p u  I  s e ) * s  i  nc ( ( 0 .  5/ p i I  * ( -  ne  u  3)  *(  pul  se- a  b  s ( t a  u  _  3 )  ) ) ; 

K3(  m,  n)  =phasel_3*phase2_3*amb_3*J  ; 

%for  receiver 

a  r  g  1  4  =(  (  Y 1 +Z 4  J .  1  )  *  (  (  u* ( ( Z4. 1 ) *  p )  )  -  (  v *  (  (  Z4 .  1  )  *  x )  )  )  ; 
phas  e 1_ 4  =e  x p( - i  *  k  y  *  a  r  g 1  4/c); 

a  r  g  2  4=(  (  Y1+Z4)  .  ,y*((p-x7+(0.  5*(u-v)*(((Yl+Z4).  1  )  *(  p+x)  )  )  ) 
p  h  a  s  e  2  4  =e  x  p ( i  *k  y  *  a  r  g  2  4/ c) ; 
tau_4=l  /  c*( ( (Y1+Z4) . 1 )*Fp-x)  )  ; 
neu_4=k  y*(  (  ( Y1+Z4) . 1 ) * ( u - v )  )  ; 
amb  4=pF*(pulse-abs(tau  4 )  )  *  e  x  p  ( - 
i  *0.  5*neu_4?pul  s  e ) *  s i  n  c ( ( 0 .  5/pi J*( - neu_4) *(  pul  se-  abs(  tau_4)  )  ) ; 

K4(  m,  n)  =phasel_4*phase2_4*amb_4!|cJ  ; 
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%  f  o  r  receiver  5 

a  r  g  1  5=(  (  Y1+Z5T.  1  ) *(  (  u*( ( Z5.  '  )*p))-(v*((Z5. ' ) *  x ) ) ) ; 
phasel  5  =e  x  p ( - i *  k  y  *  a  r  g  1  5/c); 

a  r  g  2  5=(  (Y1+Z5)  .  1  T*(  (  p  -  xj  +( 0.  5*( u- v) *( ( (Y1+Z5) . 1 )*( p+x) ) ) ) ; 
phase2_5=exp( i  *  k  _  y  *  a  r g2  5/ c)  ; 
t  au  5=1/  c*( ( ( Y1+Z5) . 1 ) *[p- x) ) ; 
neu'5=k  y * ( ( ( Y1+Z5) . 1 ) *( u- v)  )  ; 
amb_5=pf*(  pul  se-abs(tau  5 ) )  *  e  x  p  ( - 
i  *0.  5*neu_5*pul  s  e ) *  s i  n  c ( ( 0 .  5/pi J*(-neu  5)*(pul  se-abs(tau_5))); 

K5(  m,  n)  =phasel_5*phase2_5*amb_5*J  ; 

%  f  o  r  receiver  6 

a  r  g  1  6=(  (  Y1+Z6J.  1  ) *  ( (  u  * ( ( Z  6 . 1 )*p))-(v*((Z6. ' )  *x)  )  )  ; 
phas  e 1_ 6 =ex  p( - i  *  k  y  *  a  r  g  1  6/c); 

a  r  g  2  6=(  (  Y1+Z6)  .  ,T*((p-x7+(0.  5*(u-v)*(((Yl+Z6).  1  ) *( P+x) ) ) ) ; 
phase2_6=exp( i  *  k  _  y  *  a  r g2  6/ c)  ; 
t  au  6  =1  /  c * ( ( ( Y1+Z6) . 1 ) *[p- x)  )  ; 
neu~6=k  y*( ( ( Y1+Z6) . 1 ) *( u- v)  )  ; 
amb_6=pf*(  pul  se-abs(tau  6 )  )  *  e  x  p  (  - 
i  *0.  5*neu_6*pul  s  e ) *  s i  n  c ( ( 0 .  5/pi T*( - neu_6) *( pul  se-abs(tau_6))); 

K6(  m,  n)  =phasel_6*phase2_6*amb_6*J  ; 

%  f  o  r  receiver  7 

a  r  g  1  7  =(  (  Y 1 +Z  7  J  .  1  )  *  (  (  u  *  (  (  Z  7 .  1  )*p))-(v*((Z7.  1  )  *x)  )  )  ; 
phas  e 1_ 7  =ex  p( - i *  k  y  *  a  r  g  1  7/c); 

a  r  g  2  7=(  (  Y1+Z7)  .  'I*((p-xJ+(0.  5*(u-v)*(((Yl+Z7).  1  )  *(  p+x)  )  )  )  ; 
p  h  a  s  e  2  7=exp( i  * k  y  *  a  r  g  2  7/ c)  ; 
t  a  u_  7  =1/  c  *  (  (  (  Y 1  +Z7 )  .  1  )  *  [  p-  x )  )  ; 
neu~7=k  y*( ( ( Y1+Z7) . 1 ) *( u- v) ) ; 
amb  7=pF*(pul  se-abs(tau  7))*exp(- 
i  *0.  5*neu_7?pul  s  e ) *  s  i  n  c ( (  0 .  5/pi  J*( - neu_7) *( pul  se-abs(tau_7))); 

K7(  m,  n)  =phasel_7*phase2_7*amb_7*J  ; 

%  f  o  r  receiver  8 

a  r  g  1  8  =(  (  Y 1 +Z  8 1 .  1  )  *  ( (  u  *  (  (  Z  8 .  1  )*p))-(v*((Z8.  1  )  *x)  )  )  ; 
phasel  8  =e  x  p ( - i *  k  y  *  a  r  g  1  8/c); 

a  r  g  2  8=(  (  Y1+Z8)  .  'I*((p-xJ+(0.  5*(u-v)*(((Yl+Z8).  1  )  *(  p+x)  )  )  )  ; 
ph as e 2_ 8 =ex p( i  * k_ y *  a r g 2  8/c); 
t  au  8=1/  c*( ( ( Y1+Z8) . 1 ) *[p- x)  )  ; 
neu‘8=k  y*( ( ( Y1+Z8) . 1 ) *( u- v)  )  ; 
amb_8=pf*(  pul  se-abs(tau  8 )  )  *  e  x  p  (  - 
i  *  0 .  5  *  n  e  u  8*pui  sej^si  nc((0.  5/piT*(-neu  8)*(pul  se-abs(tau_8))); 

KB(  m,  n)  =phasel_8*phase2_8*amb_B*j  ; 

K(  m,  n)  =(  abs(  Kl(  m,  n)  +K2(  m,  n)  +K3(  m,  n)  +K4(  m,  n)  +K5(  m,  n)  +K6(  m,  n)  +K7(  m,  n)  +K8(  m 
,  n ) ) ) .  A2 ; 


end 

end 


%%%%%%%%%%%%%%%%% 

Localization  in  Velocity-Space 

%%%%%%%%%%%%%%%%% 


nx=l  engt  h(  u_x)  ; 
ny=l  engt  h(  u_y)  ; 
K=zeros(ny,nx); 
Kl=zeros( ny,  nx) 
K2=zeros( ny,  nx) 
K3=zeros( ny,  nx) 
K4=zeros( ny,  nx) 
K5=zeros( ny,  nx) 


%  length  of  x  =  c  o  I  u  mn  s  in  K 
%  length  of  y  =rows  in  K 
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K6  =z  e  r  o  s  ( n y ,  nx)  ; 

K7  =z  e  r  o  s  ( n y ,  nx)  ; 

K8=zeros(  n y ,  nx)  ; 

for  m= 1 :  n  y 

for  n  =1 :  n  x 

u=[u_x_start+((n-l)*res_u);u_y_end-((m-l)*res_u);0]; 

%  f  o  r  receiver  1 

a  r  g  1  1  =(  (  Y 1 +Z 1 T  .  1  )  *  ( (  u  *  ( (  Z 1 .  '  )*p))-(v*((Zl.  1  )  *x)  ) ) ; 
phasel_l=exp( - i  *k  y*argl_l/ c) ; 

a  r  g  2  1=(  (  Yl+Zl)  .  1  T*((p-xT+(0.  5*  (  u  -  v)  *  (  (  (  Y 1 +Z1 )  .  1  )  *  (  p  +x ) )  ) ) 
p  h  a  s  e  2  1  =e  x  p  (  i  * k  y  *  a  r  g  2  1/ c)  ; 
t  a  u_  1  =1/  c  M  (  (  Y 1  +Z1 )  .  1  )  *  [  p-  x )  )  ; 
n  e  u_  1  =k  y *  (  (  ( Y 1 +Z1 ) , 1 ) * ( u - v )  )  ; 
amb  l=pf*(  pul  se-abs(tau  1 ) )  * e x p ( - 
i  *0.  5*neu_l?pul  s  e ) *  s i  n  c ( ( 0 .  5/pi J*(  -  neul) *(  pul  se-abs(tau_l))); 

Kl(  m,  n)  =phasel_l*phase2_l*amb_l:,tJ  ; 

%  f  o  r  receiver  2 

a  r  g  1  2 =(  (  Y 1  +Z2J .  1  )  *  (  (  u *(  (  Z2 .  1  )* p )  )  -  (  v *(  (  Z2 .  1  )  * x )  )  )  ; 
phasel  2  =e  x  p ( - i *  k  y  *  a  r  g  1  2/c); 

a  r  g  2  2  =(  (Yl+Zl)  .  1  I  *  (  (  p  -  xl +(  0 .  5*(u-v)*(((Yl+Z2).  1  )  *  (  p  +x )  )  )  ) 
phase2_2=exp(  i  *k_y5|:arg2  2/c); 
t  au  2=1/ c*( ( ( Yl+22) .  1  ) *[p- x) )  ; 
neu~2=k  y*( ( ( Y1+Z2) . 1 ) *( u- v) ) ; 
amb_2=pr*(  pul  se-abs(tau  2 )  )  *  e  x  p ( - 
i *  0 .  5  *  n  e  u  2  * p u I  s e ) * s i  nc ( ( 0 .  5/ p i I  * ( - ne u  2)*(pul  se-abs(tau_2))); 
K2(m,  n)  =phasel_2*phase2_2*amb_2*J  ; 

%  f  o  r  receiver  3 

a  r  g  1  3=(  (  Y1+Z3J.  1  ) *(  (  u*( ( Z3.  '  )*p))-(v*((Z3. ' ) *  x ) ) ) ; 
phasel_3=exp( - i  *k  y  *  a  r  g  1  3/c); 

a  r  g  2  3  =(  (  Y 1  +Z3 )  .  1  I  *(  (  p-  xj  +  (  0 .  5  *  (  u  -  v )  *  (  (  (  Y 1 +Z3 )  .  1  )  *(  p  +x )  )  )  ) 
phase2_3=exp(  i  *k_y5|:arg2  3/c); 
t  au  3  =1  /  c * ( ( ( Yl+23) . 1 ) *[p- x) ) ; 
neu~3=k  y*( ( ( Y1+Z3) . 1 ) *( u- v)  )  ; 
amb_3=pF*(  pul  se-abs(tau  3 ) )  *  e  x  p  ( - 
i  *  0 .  5*neu_3*pul  s  e ) *  s i  nc ( ( 0 .  5/pi J*( - neu_3) *( pul  se-abs(tau_3))); 

K3(  m,  n)  =phasel_3*phase2_3*amb_3*J  ; 

%  f  o  r  receiver  4 

a  r  g  1  4  =(  (  Y 1  +Z  4  J .  1  )  *  (  (  u*  (  (  Z4.  1  )*p))-(v*((Z4.  1  )  *  x )  )  )  ; 
phas  e 1_ 4  =e  x p( - i  *  k  y  *  a  r  g  1  4/c); 

a  r  g  2  4  =(  (  Y 1  +Z4 )  ,  1  I  *(  (  p-  xj  +  (  0 .  5  *  (  u  -  v )  *  (  (  (  Y 1 +Z4 )  .  1  )  *(  p  +x )  )  )  ) 
p  h  a  s  e  2  4  =e  x  p ( i *  k  y  *  a  r  g  2  4/c); 
t  a  u_  4  =1/  c  *  (  (  (  Y 1  +Z4)  .  1  )  *  [  p-  x )  )  ; 
neu~4=k  y*( ( ( Y1+Z4) . 1 ) *( u- v)  )  ; 
amb  4=pF*(pulse-abs(tau  4 )  )  *  e  x  p  ( - 
i  *0.  5*neu_4?pul  s  e ) *  s i  n  c ( ( 0 .  5/pi T*( - neu_4) *( pul  se-abs(tau_4))); 

K4(  m,  n)  =phasel_4*phase2_4*amb_4!|cJ  ; 

%  f  o  r  receiver  5 

a  r  g  1  5  =(  (  Y 1 +Z  5  T  .  1  )  *  ( (  u  *  (  (  Z  5 .  '  )*p))-(v*((Z5.  1  )  *x)  )  )  ; 
phasel  5  =e  x  p ( - i *  k  y  *  a  r  g  1  5/c); 

a  r  g  2  5  =(  (  Y 1  +Z5 )  .  1  J  *(  (  p-  xj  +  (  0 .  5  *  (  u  -  v )  *  (  (  (  Y 1 +Z5 )  .  1  )  *(  p  +x )  )  )  ) 
phase2_5=exp( i *k_y*arg2  5/c); 
t  au  5=1/  c*( ( ( Yl+25) . 1 ) *[p- x)  )  ; 
neu‘5=k  y*( ( ( Y1+Z5) . 1 ) *( u- v)  )  ; 
amb_5=pf*( pul  se-abs(tau  5) ) *exp( - 
i  *0.  5*neu_5*pui  s  e ) *  s i  n  c ( ( 0 .  5/pi I*( - neu  5) *( pul  se-abs(tau_5))); 

K5(  m,  n)  =phasel_5*phase2_5*amb_5*j  ; 
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%  f  o  r  receiver  6 

a  r  g  1  6  =(  (  Y 1 +Z  6 1 .  1  )  *  ( (  u  *  (  (  Z  6 .  1  )*p))-(v*((Z6.  1  )  *x)  )  )  ; 
phas  e 1_ 6 =ex  p( - i  *  k  y  *  a  r  g  1  6/c); 

a  r  g  2  6  =  (  (  Y 1  +Z6 )  .  1  I  *(  (  p-  xj  +(  0 .  5  *  (  u  -  v )  *  (  (  (  Y 1 +Z6 )  ,  1  )  *(  p  +x )  )  )  )  ; 
ph  a  s  e  2_  6  =ex  p(  i  *  k_  y  *  a  r  g  2  6/c); 
t  au  6=1/  c*( ( ( Y1+Z6) . 1 ) *[p- x)  )  ; 
neu"6=k  y* ( ( ( Y1+Z6) . 1 ) *( u- v) ) ; 
amb_6=pf*(  pul  se-abs(tau  6 ) )  *  e  x  p  ( - 
i  *0.  5*neu_6*pul  s  e ) *  s i  n  c ( ( 0 .  5/pi J*( - neu_6) *( pul  se-abs(tau_6))); 

K 6 (  m,  n )  =phasel_6*phase2_6*amb_6*J  ; 

%  f  o  r  receiver  7 

a  r  g  1  7  =(  (  Y 1 +Z  7  T  .  1  )  *  (  (  u  *  (  (  Z  7 .  1  )*p))-(v*((Z7.  1  )  *x) ) ) ; 

phas  e 1_ 7  =ex  p( - i  *  k  y  *  a  r  g  1  7/c); 

a  r  g  2  7  =(  (  Y 1  +Z7 )  ,  1  I  *(  (  p-  xj  +(  0 .  5  *  (  u  -  v )  *  (  (  (  Y 1 +Z7 )  .  1  )  *(  p  +x )  )  )  )  ; 

p  h  a  s  e  2  7  =ex p( i  *  k  y  *  a  r  g  2  7/c); 

t  au  7=1/  c*( ( ( Y1+Z7) . 1 ) *[p- x)  )  ; 
neu"7=k  y*( ( ( Y1+Z7) . 1 ) * ( u - v)  )  ; 
amb  7=pf*(pul  se-abs(tau  7))*exp(- 
i  *0.  5*neu_7^pul  s  e ) *  s  i  n  c ( ( 0 .  5/pi J*( - neu_7) *( pul  se-abs(tau_7))); 

K7(  m,  n)  =phasel_7*phase2_7*amb_7*J  ; 

%  f  o  r  receiver  8 

a  r  g  1  8  =(  (  Y 1 +Z  8 1 .  1  )  *  ( (  u  *  (  (  Z  8 .  1  )*p))-(v*((Z8.  1  )  *x)  )  )  ; 
phasel  8  =ex p( - i  *  k  y  *  a  r  g  1  8/c); 

a  r  g  2  8  =(  (  Y 1  +Z8 )  .  1  J  *(  (  p-  xj  +(  0 .  5  *  (  u  -  v )  *  (  (  (  Y 1  +Z8 )  .  1  )  *(  p  +x )  )  )  )  ; 
phase2_8=exp( i *k_y*arg2  8/c); 
t  au  8=1/ c*( ( ( Y 1 +28 ) . 1 ) *[p- x)  )  ; 
neu‘8=k  y*( ( ( Y1+Z8) . 1 ) *( u- v)  )  ; 
amb_8=pf*( pul  se-abs(tau  8 ) ) *  e  x  p (  - 
i  *  0 .  5  *  n  e  u  8*pui  sej^si  nc((0.  5/piT*(-neu  8)*(pul  se-abs(tau_8))); 

KB(  m,  n)  =phasel_8*phase2_8*amb_B*J  ; 

K(  m,  n)  =(  abs(  Kl(  m,  n)  +K2(  m,  n)  +K3(  m,  n)  +K4(  m,  n)  +K5(  m,  n)  +K6(  m,  n)  +K7(  m,  n)  +K8(  m 
,  n ) ) ) .  A2 ; 


end 

end 


%%%%%%%%%%%%%%%%% 

3-D  Slice  Plot  for  Position-Space 

%%%%%%%%%%%%%%%%% 


p  x  start  20; 
p_x_end=20; 
p'y'st  art  =- 20; 
p_y_end=20; 
p'z'st  art  =- 20; 
p_z_end=20; 
r  e s _  p  =0 .  5; 

P_x=p_x_st  art :  res_p;  p  _  x  _  e  n  d ;  %  define  vector  of  x-axis  values 
P_y=P_y_st  art :  res_p;  p_  y _  e  n  d ;  %  define  vector  of  y-axis  values 
P_z-P_z_5tart :  r  e  s  _  p :  p  _  z  _  e  n  d ; 


nx=l  engt  h(  p_x)  ; 
n y  =1  engt  h(  p_y)  ; 
nz=length(p_z); 
K=zeros(  ny,  nx,  nz)  ; 
Kl=zeros(ny,nx,nz); 
K2=zeros( ny,  nx,  nz)  ; 


%  length  of  x  ^columns  in  K 
%  length  of  y  =rows  in  K 
%  length  of  z  =  pages  in  K 
%  i  n i  t  i  a  I  i  z  e  K  ma  t  r  i  x 


for  a  =1 :  n  z 

for  m=l :  n  y 
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for  n  =1 :  n  x 


p=(10A3)*[p_x  start  +(  (  n  -  l)*res_p)  ;  p_y_end  -  (  (  m- 
l)*res_p);p_z_start+((a-iy*res_p)j; 

%  f  o  r  receiver  1 

a  r  g  1  1  =(  (  Y 1  +Z 1 7 .  1  )  *  ( (  u  *  ( (  Z 1 .  1  )*p))-(v*((Zl.  1  )  *  x )  )  )  ; 
ph a s  e 1_ 1 =ex p( - i  *  k  y  *  a  r  g  1  1/c); 

a  r  g  2  1=(  (  Yl+Zl)  ,  1  T*((p-xT+(0.  5*(  u  -  v)  *(  (  (  Yl+Zl)  .  1  )  *(  p+x)  )  )  )  ; 
p  h  a  s  e  2  1 =e  x  p ( i *  k  y  *  a  r  g  2  1/c); 
t  a  u_ 1 =1/ c  *  (  (  (  Y 1 +Z1 ) .  1  )  *  [  p-  x ) ) ; 
neu" 1 =k  y*( ( ( Yl+Zl) . 1 ) *( u- v) ) ; 
amb_l=pf*(  pul  se- abs(tau_l) )*exp(- 
i  *0.  5*neu_l*pul  s  e ) *  s  i  n  c ( ( 0 .  5/pi )*(-neu_l)*( pul  se- abs( tau_l) ) ) ; 

Kl(  m,  n,  a)  =phasel_l*phase2_l*amb_l*J  ; 

%  f  o  r  receiver  2 

a  r  g  1  2  =(  (  Y 1 +Z  2  7 .  1  )  *  (  (  u  *  (  (  Z  2 .  1  )*p))-(v*((Z2.  1  )  *x)  )  )  ; 
phasel  2  =e  x  p ( - i *  k  y  *  a  r  g  1  2/c); 

a  r  g  2  2  =(  (  Yl+Zl)  .  'J*((p-xI+(0.  5*(u-v)*(((Yl+Z2).  1  )  *(  p+x)  )  )  )  ; 
phase2_2=exp( i *k_y*arg2  2/c); 
t  au  2=1/  c*( ( ( Y1+Z2) , 1 ) *[p- x)  )  ; 
neu‘2=k  y*( ( ( Y1+Z2) . 1 ) *( u- v)  )  ; 
amb_2=pr*( pul  se- abs( tau_2) )*exp( - 
i  *0.  5*neu_2*pul  s  e ) *  s  i  n  c  (  (  0 .  5/pi  j  *(  -  neu  2)  *( pul  se-abs(tau_2))); 

K2(  m,  n,  a)  =phasel_2*phase2_2*amb_2*J  ; 

K  (  m,  n,  a)  =(  abs(  Kl(  m,  n ,  a )  +  K  2  (  m,  n,a))).A2; 


end 

end 

end 


%%%%%%%%%%%%%%%%% 

Ambiguity  Implementation  for  Chirp  Signal 

%%%%%%%%%%%%%%%%% 


a  mb  _ 1 =  p  i  *(  pul  se-abs(tau_l))*exp(i  *0.  5*(  - 
neu  1  +( g  a  mma  *  t  a  u  _  1 ) )  *  p  u  I  s  e )  *  s  i  n  c  ( ( 0 .  5/pi  )  *(  -  neu_l+(  ga  mma  *tau_l))*(pulse- 
a  b  s  [  t  a  u  _  1 ) ) ) ; 

%%%%%%%%%%%%%%%%% 

Out  put  Plot 
%%%%%%%%%%%%%%%%% 


clear  all 

i  mage(p_x,  p  _  y ,  64*K/  max(  max(  K)  )  ) 
hold  on 

pi  ot(2,  0,  1  +k'  ,-5,  5,  '  Ak'  ,  15,  15,  '  Ak'  ,-2,-15,  '  Ak'  ,  15,  - 
1 5 ,  1  Ak 1  ,  1  I  i  newi  d t  h 1  ,  2 ) 
hold  off 

xl  abel ( {' \fontsi  z  e { 1 1 }  x  -  component  1  ,  1  \  fontsi  z  e  { 1 1 }  ( p  o  s i  t i  on) 1  }) 
yl  abel (  { 1 \fontsi  z  e { 1 1 }  y  -  component '  ,  '  \  fontsi  z  e { 1 1 }  ( posi  t i  on) 1  }) 
t  i  1 1  e  (  { '  \  f  o  n  t  s  i  z  e  { 1 2 }  Transmi  tter/Recei  ver  Geometry  1'}) 

%  S I  i  c  e  Plot 

nor  m=max(  ma x (  ma x (  K)  )  )  ; 

[  UX,  UY,  UZ] =meshgr i  d ( u  x,  u  y,  u  z)  ; 
si  i  c  e  (  UX,  UY,  UZ,  K/  nor  m7  [  0]  A  0] ;  [  0]  ) 
col  or  bar 


64 


LIST  OF  REFERENCES 


[1]  B.  Borden,  Radar  imaging  of  airborne  targets:  A  primer  for  applied 
Mathematicians  and  Physicists,  Taylor  &  Francis,  1999. 

[2]  M.  Cheney  and  B.  Borden,  “Imaging  moving  targets  from  scattered  waves,” 
Institute  of  Physics  Publishing  on  Inverse  Problems,  Vol.  24,  April  2008. 

[3]  K.  Tomiyasu,  “Tutorial  review  of  synthetic  aperture  radar  (SAR)  with 
applications  to  imaging  of  the  ocean  surface,”  Proceedings  of  the  IEEE,  Vol.  66, 
No.  5,  pp.  563-583,  May  1978. 

[4]  M.  J.  Gerry,  L.  C.  Potter,  I.  J.  Gupta  and  Andria  van  der  Merwe,  “A  parametric 
model  for  synthetic  aperture  radar  measurements,”  IEEE  Trans  Antennas  and 
Propagation,  Vol.  47,  No.  7,  pp.  1 179-1188,  July  1999. 

[5]  E.  G.  Larsson,  G.  Liu,  P.  Stocia  and  J.  Li,  “High-resolution  SAR  imaging  with 
angular  diversity,”  IEEE  Trans  Aerospace  and  Electronic  Systems,  Vol.  37,  No.  4, 
pp.  1359-1372,  October  2001. 

[6]  B.  Borden  and  M.  Cheney,  “Synthetic  aperture  imaging  from  high-Doppler 
resolution  measurements,”  Institute  of  Physics  Publishing  on  Inverse  Problems, 
Vol.  21,2005. 

[7]  F.  Berizzi,  E.  Dalle  Mese,  M.  Diani  and  M.  Martorella,  “High-resolution  ISAR 
imaging  of  maneuvering  targets  by  means  of  the  range  instantaneous  Doppler 
technique:  Modeling  and  performance  analysis,”  IEEE  Trans  Image  Processing, 
Vol.  10,  No.  12,  pp.  1880-1890,  December  2001. 

[8]  B.  Zheng  et  al,  “Time-frequency  approaches  to  ISAR  imaging  of  maneuvering 
targets  and  their  limitations,”  IEEE  Correspondence,  2001. 

[9]  T.  Tsao,  M.  Slamani,  P.  Varshney,  D.  Weiner  and  H.  Schwarzlander,  “Ambiguity 
function  for  bistatic  radar,”  IEEE  Correspondence,  1997. 

[10]  B.  Friedlander  and  B.  Porat,  “VSAR:  A  high  resolution  radar  system  for  detection 
of  moving  targets,”  IEE  Proc  Radar,  Sonar  Navig.,  Vol.  144,  No.  4,  pp.  205-218, 
August  1997. 

[11]  Y.  L.  Wang  et  al.,  “Space-time  adaptive  processing  for  airborne  radar  with 
various  array  orientations,”  IEE  Proc  Radar,  Sonar  Navig.,  Vol.  144,  No.  6,  pp. 
330-340,  December  1997. 


65 


[12]  B.  Borden,  “Mathematical  problems  in  radar  inverse  scattering,”  Institute  of 
Physics  Publishing  on  Inverse  Problems,  Vol.  18,  No.  1,  February  2002. 

[13]  J.  B.  Keller,  “Inverse  problems,”  The  American  Mathematical  Monthly,  Vol.  83, 
No.  2,  pp.  107-118,  February  1976. 

[14]  M.  I.  Skolnik,  Introduction  to  radar  systems,  2nd  edition,  McGraw-Hill,  2001. 

[15]  H.  D.  Griffiths,  “Developments  in  modern  synthetic  aperture  radar,”  IEEE  Radar 
Conference,  pp.  734-739,  April  2007. 

[16]  Ausherman  et  al.,  “Developments  in  radar  imaging,”  IEEE  Trans  Aerospace  and 
Electronic  Systems,  Vol.  AES-20,  No.  4,  pp.  363-400,  July  1984. 

[17]  J.  P.  Zwart,  Aircraft  recognition  from  features  extracted  from  measured  and 
simulated  radar  range  profiles,  ASCI,  The  Netherlands,  2003. 

[18]  M.  Cheney  and  B.  Borden,  Fundamentals  of  radar  imaging,  SIAM,  2009 

[19]  “Radar,”  class  notes  for  PH  4274,  Physics  Department,  Naval  Postgraduate 
School,  May  2008. 

[20]  V.  C.  Chen,  H.  Ling,  “Time-Frequency  Transfonns  for  radar  imaging  and  signal 
analysis,”  Norwood,  MA:  Artech  House,  2002. 

[21]  S.  S.  Jae,  G.  Thomas,  B.  C.  Flores,  “Range-Doppler  radar  imaging  and  motion 
compensation,”  Norwood,  MA:  Artech  House,  2001. 

[22]  M.  Bertero  and  P.  Boccacci,  Introduction  to  inverse  problems  in  imaging,  pp.  81- 
86.  Institute  of  Physics  Publishing,  UK,  1998. 

[23]  J.  R.  Guerci,  Space-time  adaptive  processing  for  radar,  Artech  House,  2003. 

[24]  J.  S.  Bergin  and  P.  M.  Techau,  “Multiresolution  signal  processing  techniques  for 
ground  moving  target  detection  using  radar,”  EURASIP  Journal  on  Applied 
Signal  Processing,  No.  47534,  2006. 

[25]  G.  Y.  Wang  et  al.,  “Dual-speed  SAR  imaging  of  moving  targets,”  IEEE  Trans 
Aerospace  and  Electronic  Systems,  Vol.  42,  No.  1,  pp.  368-379,  January  2006. 

[26]  J.  K.  Jao,  “Theory  of  synthetic  aperture  radar  imaging  of  a  moving  target,”  IEEE 
Trans  GeoScience  and  Remote  Sensing,  Vol.  39,  No.  9,  pp.  1984-1992, 
September  2001. 


66 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Professor  James  H.  Luscombe 
Code  PH/Lj 

Naval  Postgraduate  School 
Monterey,  California 

4.  Professor  Brett  Borden 
Code  PH/Bb 

Naval  Postgraduate  School 
Monterey,  California 

5.  Professor  Donald  L.  Walters 

Code  PH/We 

Naval  Postgraduate  School 
Monterey,  California 

6.  Professor  Yeo  Tat  Soon 

Director,  Temasek  Defence  Science  Institute 
National  University  of  Singapore 
Singapore 

7.  Tan  Lai  Poh  (Ms) 

Assistant  Manager,  Temasek  Defence  Science  Institute 

National  University  of  Singapore 

Singapore 

8.  Tan  Lu  Pin 
Ministry  of  Defence 
Singapore 


67 


